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PREFACE 


This  report  was  prepared  as  part  of  a  contract  effort  conducted 
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under  Task  80-13  of  Air  Force  Contract  F33615-79-C-3030.  Drs.  James  D. 
Trolingor,  William  T,  Mayo,  Jr..,  and  Anthony  E.  Smart  of  Spectron  Develop¬ 
ment  Laboratories,  Inc.,  Costa  Mesa  California,  were  the  principal 
scientists.  The  performance  period  was  April  through  November  1980.  The 
tests  were  carried  out  at  the  aerodynamic  teat  facilities  of  the  Aero¬ 
mechanics  Division  of  the  Flight  Dynamics  Laboratory,  Air  Force  Wright 
Aeronautical  Laboratories,  at  Wright-Patterson  AFB,  Ohio.  The  tests  were 
an  element  of  in-house  work  unit  24041307,  "Development  of  Testing  Tech¬ 
niques  and  Aero-Optical  Diagnostics  to  Advance  Aerodynamic  Ground  Slmula- 

'  ^  ,  '  $  v*'*  *.  k.  , 

tion"  of  Task  240413,  "Aerodypaaio1  Ground  Test  «Ta^hnplogy" . 
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SECTION  I 


INTRODUCTION 

1.  Objective 

The  objectives  of  this  study  were  several:  for  Task  I,  a  Laser 
Transit  Anemometer  (LTA)  system  was  to  be  used  for  backscatter  measure¬ 
ments  in  a  continuous  closed  circuit  tunnel,  the  Trisonic  Gasdynamics 
Facility  (TGF) .  These  LTA  measurements  were  to  be  for  characterization 
of  observable  naturally  existing  particles  existing  in  the  facility  to 
determine  particle  size  and  measurement  precision.  Under  Task  II,  veloc¬ 
ity  profiles  were  mapped  around  a  3-dimensional  microfighter  model  in 
the  TGF  and  the  Nine-Inch  Self-Adaptive  Wall  (SAW)  wind  tunnels.  The 
purpose  of  the  flow  field  mapping  program  about  the  common  rest  model  in 
both  facilities  was  for  support  of  a  program  at  Flight  Dynamics  Laboratory 
to  develop  an  adaptive  wall  concept  to  reduce  wall  interference  effects  in 
transonic  wind  tunnels.  The  9"  SAW  tun\iel  features  2D  computer  contoured 
top  and  bottom  tunnel  walls  to  correct  for  wall  interference.  The  TGF 
portion  of  the  LTA  mapping  tests  initiated  development  of  a  data  base  for 
the  caee  of  minimal  wall  interference.  The  LTA  mapping  in  the  SAW  was 
intended  to  start  evaluating  the  wall  interference  correction  when  compared 
with  the  TGF  data.  A  third  Task  (Task  III)  was  devoted  to  measurements  with 
a  thermal  emission  velocimeter  (TEV)  of  the  particle  velocities  in  the  RENT 
accelerator.  The  activity  under  Task  III,  is  reported  in  reference  1  and 
is  not  discussed  further  here.  Task  IV  of  this  program  was  the  analysis  and 
discussion  of  the  results  of  Task  I  and  II  including  recommendations  for 
future  laser  velocimeter  systems  applications  to  the  TGF  and  SAW  facilities 
and  the  requisite  data  acquisition  hardware  and  software. 
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2. 


Background 


(2  4) 

Previous  attempts  *  have  been  made  to  make  fr'n^e  laser  velo- 
cimeter  (LDV)  measurements  in  the  TGF  with  naturally  existing  particu¬ 
lates.  Observation  of  oscilloscope  traces  of  the  photomultiplier  output 
revealed  the  absence  of  particles  sufficiently  large  for  classical  burst 

(3) 

counter  detection  after  a  period  of  20  minutes  .  Attempts  to  use 

/  a  /  \ 

photon  correlation'"  ’  have  also  failed  to  produce  usable  data  after 
the  first  20  minutes  of  tunnel  run  time.  Tne  indications  have  boen  that 
the  larger  particles  (say  greater  than  .5  microns)  are  removed  by  the 
tunnel  operation  apparently  whether  or  not  the  facility  drying  and  fil¬ 
tering  systems  were  in  the  circuit.  After  the  beginning  of  tunnel  oper¬ 
ation  the  focussed  beams  are  usually  faintly  visible  aa  scatter  from  a 
large  number  of  very  small  particles  which  have  been  estimated  to  be 
less  or  equal  to  0.1  to  0.2  ym  diameter  because  of  the  visual  similarity 
of  scatter  brightness  in  the  forward  and  backwards  directions.  The 
logical  ways  co  proceed  in  such  cases  are  to  either  a)  try  to  develop 
a  reference  beam  veloclmeter  which  could  utilize  nonresolved  signals 
from  many  small  particles  at  once,  or  b)  use  a  system  such  as  an  LTA 
system  with  smaller  probe  volume  and  greater  sensitivity  to  very  small 
particles.  No  relevant  success  is  known  for  a)  above. 

Spectron  Development  Laboratories  (SDL)  has  been  developing  LTA 

(5"* 

systems  for  3  years.  Following  Schodl  '  at  the  DFVLR  in  Germany  and 
Smart ^  at  Rolls-Royce,  Smart  and  Mayo^*®^  have  conceived  and  imple¬ 
mented  at  SDL  many  detailed  improvements  of  technique  and  directed  the  pro¬ 
duction  of  several  LTA  systems  of  increasing  sensitivity,  nrecision  cf 
measurement,  and  data  processing  sophistication.  A  description  of  one 


2 


of  the  prototype  systems  fielded  in  the  summer  of  1979  is  provided  in  the 
Appendixes,  which  also  include  a  recent  (August  1980)  tutorial  discus¬ 
sion  of  LTA  principles  and  technology  with  an  extensive  bibliography. 

In  addition  to  two  versions  of  LTA  Optical  System  development  for  the 
U.S.  Air  Force  Arnold  Engineering  Development  Center,  SDL  has  fielded 
measurements  at  NASA  Ames  Research  Center,  the  General  Electric  Company, 
the  AEDC,  and  the  Naval  Postgraduate  School,  and  has  dallyered  several 
commercial  systems.  These  acti^itl^s >arfe  described  in  Appendix  A  and 
in  the  references  provided  there. 

3.  Scope 

The  following  sections  present  the  experimental  arrangements, 
the  resulting  data,  analysis  and  discussions  of  the  data,  and  recom¬ 
mendations  for  using  the  LTA  system  for  the  laser  velocimeter  system 
and  subsystems  for  the  TGF  and  SAW  tunnels.  A  summary  1b  provided 
which  refers  to  the  figures  and  other  key  results. 


SECTION  II 


LASER  TRANSIT  ANEMOMETER  USED  FOR  MEASUREMENTS 

1,  System  Description 

As  illustrated  schematically  in  Figure  1 ,  the  LTA  system  employed 
several  standard  Spectron  LTA  subsystems  and  the  recently  introduced 
Model  105  optical  head  with  150  mm  diameter  transceiver  optics.  Figure 
2  shows  the  optical  schematic  of  the  system.  Figures  3  and  A  are  photo¬ 
graphs  of  the  optical  head  and  electronics  taken  during  the  tunnel 
characterization  measurements.  A  Malvern  K7023  digital  correlator  (50 
ns  resolution)  was  used  as  the  data  processor  in  lieu  of  the  new  SDL 
Correlex  processor  (5  ns  resolution)  which  waB  not  yet  available  at  the 
time  of  the  measurements.  The  unavailability  of  the  Correlex  was  unfor¬ 
tunate  because  not  only  did  the  50  ns  correlator  have  inadequate  time 
resolution  for  meaningful  turbulence  measurements,  but  also  the  Malvern 
computer  output  circuitry  apparently  malfunctioned  producing  scrambled 
data  about  one  time  in  every  five  measurements.  The  software  of  the  LTA 
experiment  control  computer  (VGl'-lOl)  was  modified  to  ignore  occasions 
when  che  Malvern  unit  malfunctioned  and  to  repeat  the  last  measurement. 

In  Figure  1  the  MC-101  controller  and  VGU-101  videographics  unit 
include  almost  identical  subsystems  and  functional  behavior  as  the  "Data 
Management  System"  described  in  the  1979  ICIASF  paper  (see  Appendix  A) . 

The  units  used  are  simply  a  newer  version  of  the  former  "Data  Management  Sys¬ 
tem".  It  is  noted,  however,  that  the  angle  search  computation  utilized  a 
plot  of  data  rate  versus  angle  rather  than  the  data  to  data  +  background 
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Significant  Optical  C< 
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Figure  3.  SDL  Model  105  LTA  Optical  Head  Mounted  on  Traverse  at  Trisonic  Gasdynamics  Facility. 
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Figure  4.  LTA  Electronic  Data  Acquisition  Subsystem. 


figure  of  meric  described  in  the  1979  ICIASF  paper.  The  output  formats 
have  also  been  improved  as  described  below. 

The  new  SDL  Model  LTA  105  optical  head  was  utilised  in  thle  test. 
The  first  customer  for  a  Model  105  head  is  NASA  Ames  Research  Center, 
and  the  experiments  described  here  allowed  a  demonstration,  prior  to 
delivery,  of  several  extensions  of  previously  developed  technology. 

These  included  the  use  of  a  more  efficient  turning  mirror  assembly, 
more  versatile  fiber  optic  mounts,  the  new  6"  diameter,  larger  throw 
optica,  the  Lexel  Model  85  -  1.0  laeer,  improved  photomultiplier  tube 
assemblies,  and  the  Improved  RS-232  controlled  versions  of  the  discrim¬ 
inator  cards.  These  new  subsystems  all  performed  well.  The  use  of 
the  larger  optical  system  proved  to  be  advantageous  because  the  IGF 
mechanical  constraints  required  a  throw  of  30  inches.  This  much  range 
would  have  significantly  reduced  the  small  particle  sensitivity  of  the 
previously  developed  smaller  (4"  diameter)  optical  systems. 

2.  Data  Acquisition 

The  data  collection  method  for  the  experiments  consisted  of 
printing  standard  formats  1  and  2  of  the  experiment  control  software 
of  the  VGU-1Q1  (BASIC  version  used  with  Malvern  K7023  correlator). 
Notebook  logs  were  also  kept  of  position,  time,  laser  power,  and  high 
voltage  settings  for  the  photomultiplier  tube  high  voltage  power  sup¬ 
plies  (HVA  and  HVB) . 

An  example  data  prlutout  is  given  in  Figure  5.  The  top  graph 
is  a  reproduction,  with  a  peak  normalised  scale,  of  the  correlogram 
taken  with  the  angle  of  rotation  of  the  spots  set  near  the  correct  mean 
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Peak  of  Correlograra 
Normalized  to  Top  of 
Screen 


Figure  5.  Example  Data  Printout. 
11 


flow  direction.  The  bottom  graph  (on  runs  that  included  the  angle  search 
sequence)  was  the  plot  of  data  rate  vs.  angle.  A  parabolic  fit  of  the 
peak  and  one  point  either  side  of  the  peak  is  used  to  estimate  the  flow 
angle  in  situations  where  a  valid  peak  is  found  (peak  not  at  one  end) , 
the  printed  data  lists  th*  experimental  input  parameters  given  by  the 
operator  (left  hand  column)  and  the  resulting  output  values  calculated 
by  the  VGU  computer  from  the  input  parameters  and  measured  data  appear 
in  the  right  hand  column. 


Optical  System  Parameters 


The  critical  optical  system  parameters  are  listed  in  Table  1. 
Most  of  these  are  obvious  in  intent,  A  few  require  further  definition 
as  provided  here.  The  f /number  is  simply  the  standard  range/output- 
dlemeter  without  consideration  of  the  receiver  blockage  due  to  the 
turning  mirror  (see  Figure  2).  The  receiver  efficiency  is  givun  by 
the  product  of  the  (estimated)  efficiency  factors  Itemized  in  Table  2. 
The  standard  formula  for  focal  spot  diameter  would  give 


D 


(1) 


This  would  predict  a  focal  spot  diameter  of  %  12.2  pm  for  a  40  mm  diam¬ 
eter  transmitted  beam  diameter.  The  measured  value  of  16  ym  diameter 
(with  microscope  projection  technique)  reflecto  inaccuracies  in  visual 
estimates  of  and  as  well  as  any  imperfections  of  the  lenses  which 
might  tend  to  produce  larger  chan  diffraction  limited  spots.  The  laser 
power  varied  somewhat  until  it  was  discovered  thRt  an  error  in  the  wave¬ 
length  selector  prism  installation  had  been  made.  After  correction,  the 
laser  produced  300  mw  reliably. 
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TABLE  1 


OPTICAL  SYSTEM  PARAMETERS 


Description 


Symbol 


Wave /Length 

Laser  Power 

Output  Lens  Diameter 

Output  Lens  Focal  Range 

Collection  F/d  ratio 

Transmission  Efficiency 

Receiver  Efficiency* 

Tunnel  Window  Efficiency  per 
Window  per  pass 

PMT  Quantum  Efficiency 

Output  Beam  Diameter 
(1/e2  intensity) 

Focal  Beam  Diameter 

Probe  Volume  Length 

Spot  Separation 

Separation  to  Diameter  Ratio 

Detection  Filter  Impulse  Reaponse 


X 

PL 

D 

o 

R 

F/No. 


n 

w 


L 

S 

SSP. 


Value 

0.5145  ym 
100  -  300  «V 
150  mm 
745  ram 
30/6  -  F/5 
0.7 
0.4 


0.9 


0.2 


%  40  mm 
16  ym 
0.4  mm 
338  ym 
21 

25  ns 


*See  factors  listed  in  Table  2. 
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TABLE  2 


FACTORS  OF  RECEIVER  TRANSMISSION  EFFICIENCY 


1. 

Output  lens 

0.98 

2. 

Collimating  lens 

0.99 

3. 

Rotator  prism 

0.86 

4. 

Turning  mirror  and  crosa-atop 

0.65 

5. 

Pinholes 

1.00 

6. 

Fiber  optic  (diameter  loss) 

0.80 

7. 

Fiber  ends  reflection  loss 

0.92 

Total  -  (0.98) (0.99) <0.86) (0.65) (1.0) (0.8) (0.92)  -  0.40 
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The  optical  system  parameters  are  sufficient  to  determine  the 
number  of  photoelectrons  oetectable  for  a  particle  transit  in  terms 
of  the  Mie  scattering  coefficient  of  spherical  homogeneous  particles 
with  the  assumed  typical  index  of  n  ■  1.50.  This  is  done  in  the  next 
section  for  the  TGF  tunnel  characterization.  The  estimate  considers 
that  3--5  photoelectron  pulses  were  being  used  sb  the  threshold  setting. 
This  calculation  only  requires  the  detection  filter  time  constants  and 
not  the  absolute  system  gain  for  calibration. 


SECTION  II 


LASER  TRANSIT  ANEMOMETER  PRECISION  AND  SMALL  PARTICLE  SENSITIVITY 

1.  Introduction 

A  major  limitation  of  any  laser  velocimeter  working  with  naturally 
seeded  (unseeded)  air,  is  the  statistical  variability  error  which  results 
from  inadequate  measurement  precision  And  small  number  statistics.  As 
an  example  consider  the  following  which  was  typical  for  the  TGF  measure¬ 
ments:  a  single  particle  transit  produced,  for  example,  a  contribution  j 

of  one  unit  in  the  correlator  bin  number  25  of  the  48  bins  which  repre-  | 

Bent  the  correlation  function.  Since  this  is  distinguished  from  bin 

number  24  or  26,  the  single-particle  precision  in  this  example  was  4X  of  ' 

the  actual  mean  flow  (2  parts  in  48) .  An  experiment  with  zero  background  i 

light  produced  fluctuations  of  50  seconds  in  duration  with  a  mean  data  I 

pulse  rate  of  0.5/s.  In  this  example  25  transit  measurements  produce 

I 

10  units  in  store  number  25  and  15  of  store  number  26.  The  estimate  of 
delay  is  then  25.4  DT,  where  DT  -  50  ns  is  the  delay  time,  and  the 

measurement  precision  is  0.4%  as  discussed  below.  This  section  examines  < 

briefly  both  biaB  (or  systematic)  errors  and  statistical  (or  variability) 

errors  and  relates  these  to  particle  size  and  data  rate  through  system  j 

detection  sensitivity  and  probe  volume  dimensions. 

i 

2.  Bias  Errors 

Systematic  or  bias  errors  are  those  that  cannot  be  removed  by 
taking  more  data  and  averaging  long,  but  they  can  be  removed  by  better 
calibration  or  fixed  compensation..  Table  3  lists  the  investigators' 
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TABLE  3 


ESTIMATED  I.TA  SYSTEMATIC  ERROR  LEVELS 
(C.irrectahle  by  Additional  Calibration) 


Velocity  Magnitude 

1.  Spot  Separation 

2.  System  Differential  Delay 

Mach  2.3  5  ns/700  ns 
Mach  0.9  5  ns/1260  ns 

3.  Correlator  Clock  Frequency 

4.  Bias  Errors  due  to  Calculation  Algorithm 

(at  IX  flow  turbulence  or  lesa) 


Total  Accuracy  Without  Further  Calibration 
(Root  mean  square  sum) 


Flow  Angle 

Linear  Potentiometer 
Leveling  of  Optical  Head 
Parabolic  Curve  Fit 

Total  accuracy  without  filter  calibration 
(root  mean  square  sum) 


0.5!? 

0.7% 

0.4% 

<0.001% 

negligible 

1% 


0.3* 
0. 2e 
<0.1° 

0.37° 
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educated  estimates  of  the  absolute  accuracy  of  the  LTA  system  at  the 
level  of  calibration  used  in  the  tests.  These  could  all  be  improved 
by  more  careful  calibration  procedures. 


3.  Statistical  Errors  of  Velocity  Magnitude  Measurements 

In  the  recent  ISL  paper,1 a  formula  was  derived  for  the  statis¬ 
tical  error  in  the  mean  flow  magnitude  for  the  limiting  case  of  high 
detection  threshold  where  the  flat  "background"  portion  of  the  correl- 
ogram  is  zero: 


where  kQ  -  DT/tq  -  mean  (non-integer)  correlator  bin  no. 
civ/Vo  ■  Fractional  rms  velocity  deviation, 
o  /t  *  Fractional  rms  transit  delay  deviation. 

Ak  -  Fraction  of  bin  no.  that  kQ  differs  from  integer. 

N  »  Total  number  of  transit  data  pulses  crossing  both  spots. 


The  formula  thus  gives  for  the  example  of  section  3.1  with  k  -  25.4, 

o 

the  fractional  rms  error  due  to  small  number  statistics  if 


(3) 


This  is  a  result  of  0.4%  which  would  be  acceptable  in  many  applications. 

The  significance  of  a  faster  correlator  is  especially  apparent 
here.  vot  only  is  the  mean  delay  more  accurately  calculable  with  less 
interpolation,  but  the  percentage  variability  error  is  reduced  for  low 
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turbulence  flow.  With  DT  ■  5  ns  instead  of  50  ns,  then  for  the  same 
transit  delay,  kQ  *  254  instead  of  25.4.  If  the  actual  turbulence  were 
0.4%  or  less,  then  the  mean  could  be  measured  to  a  fraction  of  this 
with  only  a  few  particle  transits. 

Previous  work  had  indicated  that  setting  the  detection  threshold 
to  produce  correlogram  data  peak  to  background  ratios  of  5  to  10  seemed 
to  optimize  available  data  rate  without  introducing  excessive  additional 
variability  error  due  to  the  statistical  fluctuations  of  the  background. 
The  experimental  results  reported  in  section  IV  quickly  revealed  that 
under  conditions  of  very  little  data,  the  background-contributed 
variability  is  quite  significant  at  peak  to  background  levels  of  5  and, 
in  fact,  this  is  not  a  useful  measure  for  setting  thresholds.  All  exper¬ 
imental  measurements  after  the  Initial  teat  runs  were  made  at  low  values 
of  background  so  that  the  existing  formula  for  error  estimation  based  on 
zero  background  could  be  used. 

4.  Statistical  Errors  of  the  Flow  Angle  Estimate 

In  previous  measurements  with  low  turbulence  flows ,  the  particle 

concentration  has  been  statistically  steady  and  at  l*ast  100  data  pulse 

pairs  per  second  were  obtained  when  the  spot  pair  was  aligned  near  the 

true  flow  angle.  Under  these  conditions,  0. 1°  angle  repeatability  was 

(9) 

obtained  in  the  past.  No  detailed  theory  for  the  angular  error  has 
yet  been  developed;  however,  any  such  theory  would  have  proved  useless 
in  the  case  of  the  SAW  tunnel  resultB  due  to  particle  rate  fluctuation, 
We,  therefore,  can  only  report  experimentally  observed  variability  in 
the  angle  estimates,  which  varied  from  0.1°  to  several  tenthB  of  a 
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degree  in  the  TGF  to  more  than  one  degree  in  the  unstationary  particle 
conditions  of  the  SAW  tunnel. 


Particle  Size  Sensitivity 


This  section  presents  a  derivation  of  the  minimum  detectable  par¬ 
ticle  size  for  the  LTA  system  which  is  used  later  with  the  experimental 
results.  The  mathematical  model  used  here  is  very  simplified.  The 
statistically  expected  number  of  detected  photoelectrons  released  by 
the  photocathode  during  a  detection  integration  interval  by  a  scattering 


particle  is  <N^>  given  by 


(A) 


here  Td  is  the  transit  time  of  one  spot  (D^/V)  or  the  detection  filter 
Impulse  response  width,  whichever  is  smaller,  and  A^  is  the  rate  of 
emission  of  photoelectrons  from  the  photocathode.  If  <Nj>  is  greater 
than  3,  then  the  pulse  is  presumed  to  be  detected  by  a  discriminator 


with  the  detection  level  set  at  about  3  timeB  the  mean  photoelectron 


impulse  response  pulse  height.  Now  A^  is  calculated  as 


X 


d 


n.n  n  fti 

d  r  w  o 


(5) 


where  the  efficiency  factors  denoted  by  n  are  listed  in  Table  1,  ip  is 
the  differential  Mie  scattering  coefficient  at  180°  backscatter  (pre¬ 
sumed  constant  over  the  collecting  solid  angle  Q  for  submicron  particles 
as  an  approximation),  and  k  is  the  wave  number  2tt/A. 


vfr  K/4) 


(6) 
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and  the  collecting  solid  angle  ft  is  approximated  as 


ft 


TT 

A(F/no)2 


(7) 


combining  all  these  factors  we  obtftin: 


<V  .  Vu i zh  Av 

0  2D^  (F/no)2hv  \kz/ 

where  k  -  2n/X  -  12.21  x  10+6nfl 


(8) 

for  the  laser  radiation  wavelengths 


X  •  51A.5  nm. 


The  results  of  Equation  8  can  be  condensed  to  a  simple  result  by 
uuing  the  actual  parameter  listed  above  in  Table  1.  Thus: 

(25xlOS9) (0. 2) (0. A) (0. 7) (0. 9) A 

Nd  (2)(16xl0"6m)2(5)2(Axl0“1S(12.2xl06)2  MW  <9> 

1205  MM  • 

In  previous  work,^  we  have  modeled  the  Mie  scattering  coefficient 
for  a  lossless  spherical  scatterer  of  index  of  refraction  1.5  in  air  by  a 
simple  formula  approximately  valid  in  the  range  0.1  u®  <d£  1.25  pm,  where 
d  is  the  particle  diameter  in  meters  of 


i  (d)  -  23.7  <dxl06)\  (10) 

P 

In  order  to  determine  the  detected  number  of  photoelectrons  within 
the  impulse  width  of  the  detection  filter;  we  have  tabulated  the  results 
of  applying  Equation  (9)  with  an  assumed  laser  power  of  300  mW.  This 
tabulation  is  given  in  Table  A,  which  indicates  that  particles  between 
0.15  and  0.2  micrometers  in  diameter  were  being  observed.  It  is  pointed  out 
that  the  approximate  formula  for  the  Mie  scattering  coefficient  is  con¬ 
servative  by  a  factor  of  2  in  the  range  0.1  to  0.2  urn  diameter  which 
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TABLE  4 


PHOTON  DETECTION  NUMBER  VS.  DIAMETER 
FOR  TABLE  1  CONDITIONS 


d  (Um) 

i 

P 

N  (300mw) 
d 

N  (100  i 
d 

0.10 

0.0024 

0.86 

0.29 

0.15 

0.012 

4.32 

1.44 

0.20 

0.038 

13.7 

4.6 

0.25 

0.093 

33.5 

11.2 

0.30 

0.19 

68.4 

22.8 

0.35 

0.61 

130 

43 

0.40 

1.5 

220 

73 

0.50 

3.1 

540 

180 

23 

■  '-<•  ■  w,.- 
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allows  room  for  some  loss  on  the  part  of  the  scatterer  and/or  errors  in 
estimating  system  efficiency  factors.  For  a  vary  conservative  estimate, 
one  could  use  the  values  given  in  Table  3  for  100  mW  even  when  the  300  mW 
power  was  uaed.  The  results  still  indicate  that  0.2  yra  diameter  particles 
were  detected  at  up  to  700  m/s. 


6.  Data  Rate  and  Particle  Number  Density 

We  make  certain  deductions  concerning  the  number  of  density  of 

small  particles  from  the  LTA  system  parameters.  The  probe  volume  cross 

section  of  area  A  is  given  by: 

P 

A  -  D,L  (11) 

p  f 

where  is  the  diameter  of  each  of  the  focused  beams,  and  T,  is  the 
effective  length  of  the  probe  volume  (as  limited  by  the  receiving  optica) . 
This,  of  course,  is  only  the  true  projected  pulse-pair  capture  area  for 
particles  traveling  exactly  parallel  to  the  plane  of  the  two  focused 
beams.  For  laminar  flow  with  perfect  angular  alignment,  thBn  the  pulse 
pair  data  rate  is  given  approximately  in  terms  of  the  number  density 
of  detectable  particles  as 

R,  -  N  VA  (12) 

d  p 

which  is  easily  inverted  to  determine  number  density: 


N 


P 


Rd 

VA 


(13) 


Example  values  of  particle  number  densities  may  be  determined 
from  the  experimental  results  as  follows: 

The  Mach  0.9  data  at  about  270  m/s  gave  a  rough  average  of 
1  pulse  pair/aecond  in  the  TGF  and  50  pp/s  in  the  SAW,  These 


i 

i 
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results  correspond  to  number  densities  of  0,580  x  10^  o-' 
6  —3 

and  29  x  10  m  respectively,  of  particles  greater  than 
0.15  to  0.2  pm  diameter. 

Some  attempt  was  made  to  see  if  particle  number  density 
could  be  correlated  with  pressure,  PQ.  Although  data 
tend  to  indicate  an  increase  in  number  density  with  PQ, 

the  amount  of  data  collected  while  varying  ?  was  insuf- 

*  If  A  f  .  Ml 

ficient  to  establish  a  ojprldt  telationahip^  '  V® 

mi  w*  m  ■>  •' ,  ji .*  j 
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SECTION  IV 

EXPERIMENTAL  CHARACTERIZATION 


i 


1.  Overview 

Measurements  of  three  types  were  mede.  First  subsonic  end  super¬ 
sonic  measurements  with  and  without  seeding  were  conducted  through  the 
single  windows  in  the  2  by  2  foot  test  section  of  the  TGF  (TGF  Character¬ 
isation  Measurements) .  The  optical  setup  was  as  illustrated  previously 
in  Figure  3  with  the  system  operating  at  normal  incidence  to  the  window. 
The  second  set  of  measurements  in  the  TGF  was  mapping  the  mean  velocity 
end  flow  angle  on  selected  lines  near  a  three-dimensional  microfighter 
model  inside  the  15  x  15  inch  transonic  test  section  insert  of  the  TGF. 
This  facility  configuration  lias  two  windows  on  each  aide  as  illustrated 
in  Figure  6.  The  last  set  of  measurements  was  along  some  of  the  same 
relative  lines  near  the  same  microfighter  model  Installed  in  the  9  inch 
SAW  wind  tunnel  which  also  has  two  windows  on  oach  side,  Figure  7.  The 
mlcroflghter  data  Is  reported  in  Sections  V  and  VI.  This  section  presents 
the  results  of  the  system  characterisation  for  subsonic,  supersonic,  and 
transonic  test  suctions  with  examples  of  raw  data. 


2,  TGF  Characterisation 

The  subsonic  and  supersonic  tests  were  conducted  through  a  single 
(laminated)  Schlleren  window  2.44"  thick  with  the  LTA  beams  at  normal 
incidence  and  focused  at  a  point  approximately  at  the  center  of  the  test 
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Figure 


section  as  previously  Illustrated  in  Figure  3.  the  raw  data  consists  of 
notebook  log  entries  and  the  experiment  parameters,  correlograms ,  angle 
plots,  and  estimates  of  velocity,  sngle,  data  rate  and  background  rate 
(rate  of  building  of  correlogram  background).  Each  major  change  of  tun¬ 
nel  condition,  date,  or  purpose  was  givsn  a  new  experiment  number, 
although  run  numbers  were  incremented  throughout  all  experiments.  The 
results  are  summarized  here  according  to  experiment  number  subheading. 

a.  Experiment  1 

Wednesday,  October  22,  1980.  Shake-down  runs  with  no  data  printed 
out  by  the  Video  Graphics  Unit  (VGU) . 

b.  Experiment  2 

Mach  0.4,  stagnation  pressure  PQ  ■  2215  pounds  per  square  foot 
(PSF) .  October  23,  1980,  Runs  8-39.  During  this  set  of  data,  the  cor¬ 
relator  Interface  read  the  correlograms  with  a  shift  of  one  plaue  (each 
correlator  store  number  was  one  too  low,  so  that  the  software  interpreted 
an  actual  140.7  m/s  as  146.74  m/s,  as  for  axample  in  run  19.  For  this 
experiment,  due  to  the  low  speed,  the  discriminator  was  set  to  range  4 
(100  ns  detection  filter  gives  4  times  more  photons  within  detection 
filter  width  than  Range  2).  Laser  power  varied  between  120  and  200  mW 
during  experiment.  The  data  is  summarized  in  Table  5  for  10  consecutive 
repetitions  of  runs  at  the  best  angle.  Although  the  average  of  the 
measurements  was  within  0,63!  of  the  tunnel  computer  mean  flow  estimate, 
the  statistical  scatter  on  the  results  was  not  vary  low  (2%).  This  is 
much  more  than  would  be  predicted  by  the  theory  for  low  background  cor¬ 
relograms  with  typically  125  data  pulse  pairs  per  10  second  run.  Figure  8 
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TABLE  5 


DATA  FROM  MACH  0.4 


Natural  Seeding 
Mean  data  rate 

t‘ean  speed  (10  measurements,  10  a  each) 
RMS  Dev.  of  mean  estimates 
Tunnel  speed 
Tunnel  RMS  dev. 

P 

o 

Artificial  Seeding 
Mean  data  rate 

Run  Tunnel  Speed  m/s 

36  139.7 

37  139.8 

38  137.7 
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12.5  pulse  palra/s 
140.0  m/s 

2.73  ->  1.952 

140.8  m/s 

0.3  ■>  0.2X 

2205.0  PSF 

2200  pulse  palrs/e 
LTA  speed  m/s 

137.8 

137.8 

136.8 


illustrates  two  typical  results.  Note  that  In  the  run,  3  points  either 
side  of  the  peak  were  included  in  the  mean  flow  calculation,  and  the 
variability  thus  introduced  by  the  background  fluctuations  is  obvious  by 
inspection  of  the  two  correlograms .  Alsq  in  these  results,  the  computed 
turbulence  level  is  erratic  due  to  the  effects  of  limited  resolution. 

The  100  ns  range  of  the  correlator  was  used  due  to  the  fact  that  the  upper 
48  stores  of  the  correlator  could  not  read  out  at  all  and  had  to  be  aban¬ 
doned.  With  only  48  stores,  50  ns  mould  have  been  inappropriate.  As  an 
experiment  during  the  Mach  0.4  tunnel  condition  a  droplet  seeder  was  in¬ 
stalled  by  the  Experimental  Engineering  Branch  and  used  to  inject  a  water- 
glycerine  mixture  of  small  droplets  into  the  flow  up  stream  of  the  screens 
and  tunnel  contraction  section.  Figure  9  is  an  example  of  a  typical  result, 
and  Table  5  includes  the  corrected  new  vulocity  estimates  and  the  tunnel 
computer  estimates.  It  is  conjecture  as  to  whether  the  results  with  the 
seeder  show  systematic  error  between  the  computed  tunnel  velocity  data  and 
the  LTA  data,  or  perhaps  particle  lag  of  the  droplets.  It  is  interesting  to 
observe,  however,  that  the  correlograms  from  the  seeded  flow  indicate1 tur¬ 
bulence  levels  higher  than  one  would  suspect  for  the  tunnel  which  could  be 
due  to  particle  lag  dispersion.  The  turbulence  computed  by  the  LTA  is  not 
meaningful,  in  a  quantitative  sense,  due  to  correlator  shift  error. 

c.  Experiment  3 

Mach  0.4,  PQ  *  860  PSF.  October  23,  Runs  40-46.  The  pressure  was 

reduced  to  P  ■  860  PSF.  The  mean  data  iate  fell  to  7  per  second  (from  12.5 
o 

at  Pq  »  2205  PSF).  A  better  procedure  for  optimum  data  accuracy  might  have 
been  to  increase  the  detection  threshold  to  reduce  the  relative  amount  of  the 
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2SP0T  SEP<MM> 
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14VE 
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NOTE:  Correlator  shift  error  not  included  here. 


Figure  9.  Mach  0.4  ReeultB  with  Droplet  Seeding. 
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background  to  signal.  The  results  of  3  ten  second  measurements  yielded  an 
average  estimate  of  139.8  m/s  with  a  standard  deviation  of  3  m/s.  The  aver¬ 
age  of  two  100  second  measurements  was  136.0  m/s  with  a  difference  between 
the  two  of  1.6  m/s.  The  tunnel  computer  indicated  a  mean  speed  of  138.7  m/s 
with  a  standard  deviation  of  0.73  m/s  as  determined  from  conventional  pressure 
and  temperature  measurementa .  The  agreement  is  actually  remarkable  consider¬ 
ing  the  very  low  quality  of  the  correlograms  due  to  Incorrect  threshold  setting 
illustrated  in  the  right  side  correlogram  in  Figure  10. 

d.  Experiment  4 

Mach  0.7,  P  ■  1319  PSF.  Runs  47-59,  October  23,  1980.  The  previous 
o 

results  had  indicated  the  need  for  a  change  of  threshold.  This  is  a  relative 
setting  which  can  also  be  effected  by  changing  the  PMT  high  voltage  (HV) 
settings.  Decreasing  the  PMT  HV  Is  somewhat  equivalent  to  increasing  the 
detection  threshold  within  limits.  This  is  limited  by  the  constant  2.2 
volt  overload  detection  threshold  which  rejects  large  signals  by  inhibiting 
the  discriminator  output,  and  on  the  other  extreme,  by  drift  and  noise  of 
the  discriminators  at  the  level  of  a  few  millivolts.  The  detection  threshold 
range  is  generally  maintained  at  100  to  200  mv.  Generally  speaking,  the 
pulse  detection  rates  on  discriminator  channels  A  and  B  were  kept  matched 
by  using  the  correlator  monitor  stores  to  check  balance  when  adjusting  the 
HVA  and  HVB  settings.  A  fault  in  the  Malvern  computer  output  circuit  caused 
Incorrect  readings  at  times  to  be  printed  in  the  data  for  Rate  A  and  Rate  £ 
so  that  the  apparent  disparity  of  the  printed  data  is  not  factual. 

By  lowering  the  PMT  gain  and  thus  effecting  a  relative  increase 
in  the  detection  thresholds,  the  typical  correlograms  for  the  Mach  0.7 
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runs  were  made  to  appear  as  shown  in  Figure  11,  where  the  right  aide 
correlogram  illustrates  the  effects  of  a  higher  threshold  setting. 

Table  6  summarises  the  results  of  S  consecutive  rune  of  5  seconds 
interval  and  also  two  runs  of  50  seconds  each.  Even  though  3  points 
either  side  of  the  peak  of  the  corrtlogram  were  still  used,  the  measure¬ 
ment  standard  deviation  was  reduced  by  the  increase  in  the  relative 
detection  threshold.  Also,  it  was  observed  that  the  variation  in  the 
total  number  of  data  pulse  pairs  per  run  significantly  exceeded  the 
predictions  of  Poisson  statistics,  as  evidenced  by  the  sequence  of  data 
rates  given  in  the  table. 

It  was  observed  that  the  average  of  the  5  consecutive  5  second 
runs  agreed  with  the  tunnel  estimate  average  absolutely  to  within  0.3X, 
This  may  be  a  fortunate  accident  since  the  average  of  the  two  50  second 
runs  only  agreed  to  within  0.8X.  However,  for  absolute  agreement  with¬ 
out  collaboration  or  extensive  calibration,  this  seems  quite  good. 

e.  Experiment  5 

Mach  2.3  PD  •  715  PSF.  Runs  60-64.  October  23,  1980.  This 
series  was  started  with  only  25  minutes  until  tunnel  shut-down  time. 
There  were  difficulties  due  to  the  Malvern  readout.  With  limited  test 
time  remaining  for  the  day  it  was  decided  to  use  the  seeder  to  Increase 
data  rate  for  remaining  runs.  The  results  are  shown  in  Figure  12.  The 
velocities  wore  corrected  by  the  one  store  error  as  shown  in  the  figure. 
The  repeatability  was  to  0.12X.  The  tunnel  computer  gave  545.7  m/s, 
which  differ*  by  3X.  The  source  of  the  discrepancy  i#  not  known,  but 
until  shown  otherwise  both  tunnel  calibration  and  the  LTA  system  are 
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TABLE  6 


Run 

Data 

MACH  0.7 

Po  -  1319  PSF  DATA  SUMMARY 

Duration 

Race 

(s'1) 

Corrected 

Velocity 

m/a 

Tunnel 

Estimate 

m/a 

48 

6 

234.2 

236.7 

5s 

49 

3 

234.0 

236.8 

5b 

50 

16 

235.1 

237.1 

5b 

51 

6 

239.5 

N/A 

5b 

52 

13 

237.2 

236.1 

Mean, 

Dev. 

8.8 

.5.4 

236.0,2.3 

236.7,0.42 

57 

5 

235.3 

50s 

58 

4 

-  _  - 

234.4 

236.7 

30b 

Mean , 

Dlff . 

4.5 

.1 

234.8,0.9 

*V  -  563  1 
Duration  m 


5  s 


See  Figure  5  For  Detailed  Explanation. 


Figure  12.  Mach  2.3  Data  with  Seeding. 


suspect.  Although  there  should  have  been  no  bias  errors  due  to  the  some¬ 
what  coarse  correlator  time  resolution,  it  would  certainly  have  been  pre¬ 
ferable  to  use  the  5  ns  correlator  for  the  measurements. 

( .  Experiment  6 

Mach  2.3,  Pq  -  2846  PSF,  Friday,  October  24,  1980,  Runs  65-88.  A 
series  of  runs  was  made  while  the  tunnel  was  approaching  condition  and 
after  reaching  condition.  The  data  rate  varied  widely  as  did  the  measure¬ 
ments,  as  illustrated  in  Table  7  and  Figure  13  where  the  left  hand  side 
represents  the  low  data  rate  case  and  the  right  hand  side  represents  the 
high  data  rate.  Reasonable  agreement  with  the  tunnel  estimate  of  560.1  m/s 
was  obtained  when  the  data  rate  was  adequate,  but,  it  is  clear  from  Fig¬ 
ure  13  and  Table  6  that  several  things  could  have  been  improved.  The 
threshold  was  still  sufficiently  low  to  allow  significant  background  when 
the  data  rate  was  low.  This,  combined  with  inadequate  correlator  time 
resolution  and  the  variation  of  2  points  of  the  correlogram  either  side 
of  the  peak, resulted  in  lack  of  repeatability.  Near  the  unplanned  shut¬ 
down  of  this  series  (due  to  a  tunnel  cooling  system  problem)  another 
attempt  at  seeding  was  made.  Runs  85  and  86  did  not  show  the  same  large 
data  rates  obtained  the  day  before,  presumably  because  the  seeding  stream¬ 
line  missed  the  probe  volume. 

Just  prior  to  tunnel  shut-down  it  was  observed  that  the  tunnel  con¬ 
tained  a  very  fine  mist,  clearly  visible  by  scattering  from  the  laser 
beams,  however,  since  there  was  no  large  data  rata,  these  droplets  were 
smaller  than  the  detection  size  (estimated  at  £  0.2  pm  diameter). 
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The  axial  flow  direction  is  defined  to  be  270  degrees  with  the  convention  that  0  degrees 
is  vertical  and  above  the  centerline  of  the  tunnel. 
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g.  Experiment  7 


Microfighter  Transonic  Mach  0.7,  Pq  m  3218  PSF,  October  27,  1980, 

Runs  83-102.  The  micrcfighter  model  was  installed  inverted  in  the  TGF 
15-inch  transonic  teat  section.  The  setup  is  illustrated  in  Figure  6. 

The  outer  window  was  the  same,  2.44"  thick,  and  the  inner  window  1.22" 
thick  with  3.27"  separation  between  the  windows.  The  spot  separation 
was  remeasured  by  the  fringe  method  with  both  windows  In  place,  A  light 
change  is  within  the  measurement  precision  of  the  spot  separation. 

The  experiment  control  software  was  revised  to  correct  for  the 
displacement  by  one  store  which  had  occurred  in  the  previous  measure¬ 
ments.  An  interim  review  of  the  subsonic  and  supersonic  data  revealed 
the  need  for  higher  thresholds,  even  at  the  expense  of  lost  data  when 
data  rates  are  low.  Also,  additional  software  changes  were  made  to  test 
each  correlogram  for  consistency  of  numerical  levels  of  the  background 
with  previous  data  and  rerun  any  measurements  for  which  the  Malvern  unit 
produced  unusable  data.  This  allowed  successful  angular  measurements  to 
be  made  (during  the  previous  measurements  the  probability  of  obtaining  5 
successive  usable  correlograms  from  the  correlator  was  too  small  to  be 
useful) . 

In  order  to  characterise  the  precision  obtained  experimentally,  a 
series  of  rune  was  made  at  a  location  4"  below  the  model  on  the  upstream 
end  (near  free  stream  conditions).  The  mean  velocity  and  direction  were 
repeatedly  measured  from  the  time  the  tunnel  began  running.  Due  to  a  higher 
scattered  or  flare  light  levels  from  the  two  windows,  the  PMT  high  voltage 
had  to  be  turned  down  about  100  volts  with  respect  to  earlier  measurements 
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in  order  to  obtain  low  background  correlograms.  This  resulted  In  a 
lower  overall  data  rate. 

Table  8  is  a  summary  of  the  computed  results  of  this  series. 

The  measured  velocity  (bottom  graph)  and  angle  history  (top  graph)  are 
also  plotted  In  Figure  14.  The  mean  and  standard  deviation  of  the 
measured  magnitude  as  compared  with  the  tunnel  estimates  were  ao  follows: 


Mean 

S,  D. 

S.  D.% 

LTA 

233.1  m/s 

2.67 

1.15% 

Tunnel  Computer 

236.45  m/s 

0.83 

0.35% 

Angle  LTA 

270.40° 

0.31° 

These  results  are  actually  quite  amazing  when  one  considers  that 
at  the  lowest  data  rate  about  1  per  second ,  the  correlograms  for  a  5 
second  run  have  about  5  correlated  pulse  pairs  in  a  background  level  of 
less  than  one  per  store  as  illustrated  by  the  left  hand  side  of  Figure  IS. 
Run  119  shown  on  the  right  hand  side  of  Figure  IS  was  taken  with  50  seconds 
per  angle  to  observe  the  effects  of  more  deta.  The  result  is  a  measurement 
much  closer  to  the  average  measurement  than  tha  standard  deviation  of  the  5 
second  data  sets.  We  note  that  the  LTA  measurement  is  a  little  below  that 
of  tha  tunnel  computer.  This  is  probably  the  result  of  a  combination  of 
both  statistical  fluctuations,  the  difference  between  any  systematic  errors 
of  both  the  tunnel  calibration  and  the  LTA  and,  of  course,  any  actual 
reduction  in  velocity  at  the  au'ssur ament  location  due  to  the  presence  of 
the  model. 

At  the  end  of  the  sequence  of  runs  just  described,  the  LTA  probe 
volume  was  moved  up  near  the  nose  of  the  model  in  an  attempt  to  find  the 
stream  line  of  droplets  produced  by  the  seeder.  The  data  rate  increased 
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TABLE  8 

SUMMARY  OF  COMPUTED  DATA 


Run 

No. 

Tine 

S 

Data 

Rate 

s'1 

Recorded 

Velocity 

M/S 

Angle 

Deg. 

Duration 

Total 

Data 

Pulse  Pairs 

93 

11:36:37 

10 

403.6 

270.45 

58 

50 

P 

94 

11:40:32 

36 

418.3 

271.57 

5s 

180 

P 

95  • 

11:44:36 

6 

425.3 

270.9 

5s 

30 

P 

96 

11:47:56 

3 

417.9 

271.64 

58 

15 

P 

97 

11:50:56 

18 

260.6 

270.83 

5b 

90 

P 

J 

98 

11:55:09 

2 

235.3 

270.90 

5s 

10 

P 

99 

11:58:47 

7 

235.3 

270.90 

58 

35 

P 

100 

12:02:00 

3 

231.1 

269.78 

5s 

15 

P  Warm 

101 

12:05:31 

6 

235.9 

270.83 

5b 

30 

P  Up 

102 

12:11:10 

5 

236.6 

270.45 

5a 

25 

P 

on 

condition 

103 

12:15:25 

4 

229.0 

270.16 

5a 

20 

P 

104 

12:19:10 

3 

230.1 

270.44 

5b 

15 

l> 

107 

12:25:00 

2 

231.6 

270.39 

5b 

10 

P 

109 

12:29:30 

2 

235.6 

270.36 

5b 

10 

P 

110 

12:31:45 

3 

228.9 

270.54 

5s 

10 

P 

112 

12:35:39 

1 

236.2 

270.52 

5* 

5 

P 

113 

12:38:23 

2 

232.2 

270.45 

5b 

10 

P 

115 

12:44:00 

1 

229.1 

270.47 

5b 

5 

P 

116 

12:49:10 

1 

234.7 

269.89 

5s 

5 

P 

119 

1:01:33 

1 

232.1 

270.67 

50s 

50 

P 

120 

1:35:05 

45 

228.1 

283.7 

5s 

225 

S 

121 

1:37:44 

8 

228.4 

283.78 

5b 

40 

P  s 

122 

1:40:51 

7 

226.9 
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Figure  15.  Very  Low  Data  and  Background  Correlograns  from  Tr?nsnrii-  Flows. 


co  45/second,  but  was  nowhere  near  Chat  obtained  in  the  earlier  singes 
window  measurements  with  Che  seeder.  This  could  be  due  to  the  fact  that 
the  PMT  high-voltage  was  reduced  due  to  increased  flare  light  or  due  to 
poor  streamline  probe  volume  relative  positioning.  After  some  trials, 
it  was  decided  that  the  seeding  streamline  could  not  be  readily  moved 
far  enough  down  to  reach  the  location  4"  below  the  model  centerline  at 
which  the  primary  measurements  were  desired. 

The  tunnel  was  shut-down  and  several  changes  were  made.  The  seed¬ 
ing  probe  was  modified  in  an  attempt  to  allow  seeding  et  4"  below  the 
model  centerline,  and  the  ITA  system  waa  inclined  at  angle  of  6s  to  the 
tunnel  window  (rotation  about  a  vertical  axis)  in  order  to  reduce  window 
flare  light.  This  last  measure  proved  helpful  and  showed  that  there  were 
multiple  reflection  flare  effects  produced  by  going  through  two  seta  of 
windows  which  were  not  sufficiently  rejected  in  some  positions.  Anti¬ 
reflection  coating  on  the  inner  window  or  a  movable  flare  stop  between 
the  window  would  have  been  helpful.  By  using  a  6*  angle  of  incidence, 
all  of  the  multiple  window  effects  were  eliminated  except  light  from  the 
inner  window  which  was  reflected  and  focused  on  the  outer  window.  This 
problem  occurred  at  only  a  few  very  specific  probe  volume  depths  and  did 
not  affect  most  of  the  succeeding  measurements. 

h.  Experiment  8 

Microfighter  Transonic  Mach  0.7,  ■  3230  PSF,  October  28,  1980, 

Huns  124-149.  After  the  modifications  were  mode,  a  series  of  runs  were 
conducted.  Despite  the  improvements  which  allowed  raising  the  PMT  high 
voltage  levels  to  even  higher  than  was  used  with  the  single  windows,  the 
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data  pula e-palr  rate  remained  at  one  per  second  or  lesa.  This  fact 
required  the  measurement  duration  to  be  set  at  50  seconds.  Measure¬ 
ments  were  made  5  successive  times  at  a  point  4"  below  the  model 
centerline  (runs  131-135),  2  times  at  a  point  1/16"  above  the  model 
centerline  just  above  the  surface  of  tho  model,  and  again  4  measurements 
2"  below  tha  centerline  of  the  model.  The  results  in  Table  9  may  be  re¬ 
lated  to  the  coordinates  x,y,s  given  in  the  next  section. 

These  results  wsre  the  last  of  the  characterisation  measurements 
before  the  microfighter  flow  mapping  was  undertaken. 

A  few  additional  modifications  were  made  to  the  equipment  on 
Wednesday,  October  29.  These  consisted  of  changing  the  6*  off-axis 
angle  of  the  instrument  with  respect  to  the  window  to  that  of  pitch  (down) 
instead  of  yaw  to  sea  if  data  rate  was  affected  (it  was  not).  Also,  opti¬ 
cal  stop  was  inserted  in  the  LTA  heed  to  eliminate  one  additional  lena 
flare  source.  However,  these  changes  had  little  affect  on  the  system 
performance. 
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TABLE  9 


MEASUREMENTS  NEAR  THE  MICROFIGHTER  AT  Z»0.  X— 1.5 


Y 

Maan  V 

M/S 

SDV% 

Maan 

Angla 

Dag. 

Angla 

SD 

-4.0 

251.5 

0.5% 

272.3 

N/A 

+11/16 

255.5 

0.3% 

258.2 

N/A 

-2 

235.0 

0.9% 

270.4 

0.23* 

51 


SECTION  V 


AERODYNAMIC  MEASUREMENTS  IN  THE  TRISQNIC  GASDYNAMICS  FACILITY 

1.  Overview 

Figures  16a  and  b  are  three-fourths  scale  drawings  of  the  micro¬ 
fighter  model  ahowing  the  model  coordinate  system  used  aa  a  reference  system 
for  the  TGF  and  SAW  wind  tunnel  testa.  The  model  was  mounted  at  eero  angle 
of  attack  for  the  entire  teat  aeries  in  both  tunnels.  It  was  mounted  upside 
down  in  the  TGF  and  right  side  up  in  the  SAW  wind  tunnel.  Data  were  acquired 
on  0.5"  increments  along  lines  located  above  and  below  the  model  centerline 
at  y  ■  ±4.0".  The  precise  location  of  the  chosen  origins  and  axes  is  pre¬ 
sented  in  subsequent  figures  with  representative  data. 

To  check  consistency,  some  of  the  data  pointB  ware  repeated  on  differ¬ 
ent  days,  as  Illustrated  on  data  plots.  The  dat< rates  obtained  were  gener¬ 
ally  from  about  0.5  to  1  per  second.  For  these  low  data  rates,  50  seconds 
per  spot  pair  angle  were  used  to  assure  good  quality  correlograms.  The 
experlmentel  run  numbers  wore  150-270  which  comprise  experiment  numbers  9-15. 
Figure  17  illustrates  typical  results  after  the  laser  power  was  adjusted  to 
300  mW.  The  angle  plot  shown  in  Figure  17  is  typical.  In  some  cases  the 
statistical  fluctuations  of  data  rate  with  time  are  clearly  limiting  the 
angle  estimate  precision.  When  the  data  output  format  software  was  written, 
it  was  not  envisioned  that  rounding  the  data  and  background  rates  to  whole 
units  would  pose  any  limitations.  The  data  to  background  ratio  does  retain 
additional  decimal  places  which  provide  a  more  accurate  measure  of  the  back¬ 
ground  rate. 


■ 1  ■  *•  it 
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NOTEj  See  Figure  5  for  detailed  explanation. 


Figure  17.  Example  of  Data  from  TGF  Microfighter  Test. 
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2. 


Measurements  at  Mach  0.7  and  0.9 


Figure  IS  shows  plots  of  mean  velocity  versus  position  along 
lines  at  4  Inches  above  and  below  the  model  at  two  Mach  numbers.  Figure  20 
contains  the  corresponding  plots  of  flow  angle.  Scans  across  the  tunnel  at 
y  ■  -4,  x  •  0  were  also  made.  Results  of  these  scans  are  presented  In  Fig¬ 
ures  21  and  22. 

It  was  observed  that  the  tunnel  computer  printout  of  velocities  during 
these  Mach  0.7  runs  had  a  mean  value  of  237.57  m/s  and  a  standard  deviation  of 
0.62  m/s  over  the  35  runs  of  Experiment  9.  These  tunnel  values  are  Indicated 
on  Figure  18  at  the  left  hand  side  of  the  figure. 

In  order  to  experimentally  observe  the  repeatability  of  the  measure¬ 
ments,  portions  of  the  centerline  scan  at  y  -  -4"  were  repeated  on  different 
days  as  Indicated  in  the  figures.  The  small  group  of  5  measurements  at  x  -  0 
were  taken  near  the  beginning  of  tunnel  operations  when  the  data  rate  was 
higher  than  average;  hence  the  smaller  statistical  scatter  in  these  points. 

The  mean  and  standard  deviation  of  the  tunnel  computer  printout*  for  the  free 
stream  velocity  are  shown  on  Figure  18  for  comparison. 

3 .  Discussion  of  Results  in  the  TGF 

The  mean  velocity  measurements  in  the  TGF  generally  had  statistical 
scatter  of  IX  or  less.  These  results  were  obtained  typically  with  10-50 
data  pulse  pairs  obtained  in  a  50  second  data  collection  interval.  The 
results  indicate  that  a  faster  correlator  wnuld  have  been  helpful  as  would 
the  availability  of  even  a  small  amount  of  additional  scattering  particles. 
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Figure  13.  Velocity  Measurements  at  Mach  0.7  and  0.9  (Continued) 
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The  turbulence  estimation  algorithms  which  are  included  in  the 
LTA  software  package  were  not  useful  in  this  exercise  due  to  the  coarse 
correlator  resolution  and  the  very  small  amounts  of  data.  A  faster 
correlator  would  presumably  have  allowed  such  measurements. 


SECTION  VI 


MICRCFIGHTER  MEASUREMENTS  IN  THE  SELF-ADAPTIVE  WAU  TUNNEL 

1.  Tunnel  Characterization 

After  the  completion  of  the  TGF  measurements ,  the  LTA  system  was 
moved  to  the  Self-Adaptive  Wall  (SAW)  wind  tunnel  which  has  a  9  Inch 
square  test  section  and  operates  in  a  blowdown  mode.  It  was  desired  that 
the  same  array  of  data  points  be  measured  in  the  SAW  aB  had  been  measured 
in  the  TGF  for  the  purpose  indicated  In  the  Introduction  (page  1) .  How¬ 
ever,  it  should  be  noted  that  the  stagnation  pressure  was  35  PSIA  for  both 
Mach  numbers  0.7  and  0.9.  This  was  substantially  higher  than  in  the  TGF. 

In  the  SAW  wind  tunnel  the  model  was  to  be  oriented  in  the  more  usual 
fashion  (i.e.,  cockpit  up),  so  the  supports  for  the  LTA  Bystem  were  changed 
from  6°  to  4°  to  avoid  incidence  of  the  laser  beams  on  the  translucent  white 
silicone  rubber  rods  which  comprise  the  top  walls?  the  prime  region  of 
interest  of  the  measurements  was  y  -  -4  and  z  ■  0.  Figure  22  is  a  schlieren 
photograph  of  the  microfighter  taken  in  the  nine  inch  SAW  wind  tunnel. 

Prior  to  commencing  a  series  of  runs  to  duplicate  the  TGF  points,  a 
few  runs  were  made  to  check  out  the  system  and  the  SAW  tunnel  character¬ 
istics.  No  Internal  changes  or  alignments  were  required  of  the  LTA  optical 
head  after  the  move  from  the  TGF. 


*At  the  limits  of  sensitivity  light  colored  walls  which  partially  diffuse 
and  retransmit  light  have  been  found  to  be  detrimental  to  sensitivity 
due  to  the  increased  background  photon  level.  Opaque  low  reflectance 
surfaces  are  preferable. 


j„it-  .1  -  .  _ , _ r-.-—  uAil  ••  ■>» WW***-1*  wMJ* 

jf  MKCUttNO  PX<M  BUMUMOT  FILMED 
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Figure  22.  Schlieren  Photo  of  the  Microfighter  Model  in  the  SAW  Tunnel. 


These  first  checkout  SAW  runs  produced  some  interesting  results. 

The  data  rate  was  as  high  as  50  per  second  on  some  5  second  runs.  An 
example  of  a  typical  printout  is  shown  in  Figure  23.  This  was  excellent 
because  it  meant  that  100  to  250  data  pairs  could  be  obtained  in  a  5 
second  run.  As  illustrated  earlier  in  Section  III  ,  100  data  pulse  pairs 
would  typically  give  rms  deviation  of  leas  than  0.3X.  The  results  in 
the  l'GF  measurements  typically  exhibited  scatter  of  leas  than  IX  with 
fewer  than  25  data  pulse  pairs.  This  is  further  evidence  that  repeatabil¬ 
ity  of  measurements  with  no  more  than  0.5%  rms  deviation  would  be  expected 
in  the  SAW. 

The  next  interesting  observation  concerned  flare  light.  It  was 
observed  that  when  a  reflection  from  the  inner  window  of  the  transmitted 
beam  comes  to  focus  on  one  of  the  leas  than  perfectly  clean  surfaces  of 
the  outer  window,  light  scattered  by  the  outer  window  surface  will  reflect 
back  (by  way  of  the  inner  window  surface)  on  to  the  LTA  photo  detectors. 
There  is  no  way  to  shield  or  protect  the  photo  detectors  from  this  flare 
light  source  from  inside  the  LTA  head.  The  problem  can  be  minimized 
externally  aB  follows: 

a)  Clean  all  surfaces  of  outer  windows  carefully  to  reduce 
scatter, 

b)  If  feasible,  antireflection  coat  the  inner  windows, 

c)  If  possible,  direct  the  input  beam  at  an  angle  to  the 
window  and  use  a  reflection  atop  not  located  at  the  focus 
of  the  reflected  transmitted  beams,  or 

d)  Simply  avoid  the  isolated  locations  for  which  the  reflection 
focus  occurs. 
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Figure  23.  Typical  5  Second  Data  Run  From  SAW  Tunnel. 
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1'he  third  item  of  interest  observed  on  the  checkout  runs  was  the 
highly  time  variable  scattering  particle  concentration.  The  oscillo¬ 
scope  traces  of  the  raw  signal  pulses  showed  visually  changing  pulse  pair 
data  rates  which  varied  noticeably  at  random  intervals  from  a  fraction  of 
a  second  to  several  seconds.  Next,  visual  observation  of  the  laser  beams 
during  some  of  these  checkout  runs  showed  that  the  light  scattered  from 
the  beams  as  a  whole  varied  greatly,  sometimes  in  "puffB"  end  sometimes 
more  slowly  (seconds'  duration).  These  observations  ware  made  during 
several  runs,  each  of  from  1.5  to  1.75  minutes  duration,  in  which  5  to  7 
different  measurements  of  5  seconds  each  were  made  in  an  angle  search 
sequence.  Because  the  angle  estimation  algorithm  in  its  present  form  is 
based  on  the  assumption  of  stationary  particle  concentration  statistics, 

these  findings  ruled  out  the  use  of  the  instrument  for  flow  angle  measure- 

A 

ments  without  modifications  of  the  angle  algorithm.  No  records  were  obtained 
of  the  worst  examples  of  the  effects  of  particle  rate  fluctuation  due  to  the 
fact  that  the  "peak"  was  found  at  the  extreme  end  of  the  angle  sequence, 
which  the  software  considers  as  "no  peak  found"  and  did  not  record.  How¬ 
ever,  Figure  24  is  one  example  of  the  results  of  an  angle  search. 


*  The  moat  sensitive  way  to  use  the  LTA  system  is  to  set  the  detection 
thresholds  low  enough  to  detect  many  single  photo-electron  pulses  in 
the  presence  of  only  a  few  correlated  pulse  pairs  and  evan  when  the 
threshold  was  "high",  several  times  as  many  background  photon  pulses 
as  signal  pulse  pairs  were  detected.  When  data  rata  and  small  particle 
sice  are  not  a  problem,  then  by  turning  the  detection  threshold  suffi¬ 
ciently  high,  only  particle  produced  pulses  will  be  detected.  Under 
such  conditions,  the  number  of  data  versus  angle  can  be  normalized  by 
the  number  of  pulae#  on  channel  A  to  obtain  the  conditional  probability 
that  a  particle  ia  detected  by  channel  B  if  it  is  detected  by  channel  A. 
This  is  nothing  new.  It  has  been  used  by  both  Smart  at  Rolls  Royce  and 
Schodl  at  DFVLR  in  the  peat.  But  it  does  not  allow  the  full  sensitivity 
of  the  instrument  to  be  utilized,  and  this  approach  has  not  been  uBed 
recently. 
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Figure  24.  Example  of  an  Angle  Search  In  the  SAW  Tunnel: 
Data  VerauB  Angle  with  5  Seconds /Angle. 
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2. 


Measurements  Near  the  Microfighter  Model 


After  the  checkout  runs  additional  tasting  was  devoted  to  making  a 
series  of  measurements  along  two  lines  above  the  microfighter  model  at 
Mach  0.9  for  velocity  only.  Flow  angle  measurements  could  not  be  made 
for  reasons  described  on  the  previous  cage.  The  same  spot  orientation 
angle  was  used  throughout  (270  degrees  is  parallel  to  the  tunnel  center- 
line).  This  is  of  little  consequence  to  the  mean  flow  velocity  magnitude 
measured  in  this  manner  because  there  is  high  probability  the  maximum  flow 
angle  deviation  is  no  more  than  2*  from  horizontal  and  the  cosine  of  2s 
is  0.9994. 

Some  results  of  the  SAW  measurements  are  plotted  in  Figure  25.  Several 
repeat  and  fill-in  points  were  taken  in  an  attempt  to  determine  the  repeat¬ 
ability  (statistical  fluctuation  level)  of  these  measurements.  It  became 
apparent  that  the  tunnel  velocity  conditions  were  difficult  to  maintain  and 
repeat  to  the  precision  already  established  for  the  LTA. 
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Figure  25.  Mean  Velocity  Along  the  Line  Z 


SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

Aft«r  running  the  TGF  for  20  to  30  minutes^  tha  flow  contalna  on 
the  average  laaa  than  0.5  particle  per  cc  greater  than  0.2  pm  In  diameter. 
These  teata  still  have  not  defined  the  number  which  exist  below  0.2  ym 
in  diameter  in  question.  Thia  produces  1  or  leas  data  pulse  pairs  per 
second  at  transonic  speeds  rith  the  SDL  LTA  system.  Previous  experience 
shows  that  fringe  LV  systems  with  larger  probe  volumes,  which  would  inter¬ 
cept  more  particles  per  second,  have  not  had  adequate  sensitivity  to  be 
usable.  One  conclusion  is  that  there  are  too  few  of  the  larger  particles 
(nearer  to  1  pin  in  diameter)  that  would  be  needed  for  a  fringe  system  to 
operate  after  initial  start  up. 

The  LTA  system  was  demonstrated  to  be  capable  of  making  accurate 
backscatter  mean  flow  speed  and  angle  measurements  in  a  2  foot  test 
section  of  the  TGF  without  seeding  at  all  Mach  numbers.  In  a  transonic 
test  section  measurements  were  made  through  2  windows  without  artificial 
seeding  at  a  rate  of  1  measurement  per  6  minutes  and  at  an  error  level 
of  less  than  IX.  The  results  indicate  that  a  faster  (5  ns)  correlator 
and  faster  data  storage  (without  printing  each  run)  could  have  obtained 
error  levels  of  a  small  fraction  of  one  percent  in  mean  flow  in  less 
than  two  minutes  per  point  (by  reducing  the  data  collection  time  from  50 
seconds/angle  to  about  10  to  15  seconds). 

Since  t>-  .imally  designed  LTA  system  can  respond  to  at  least  0.5 
micron  particles,  it  should  be  a  simple  matter  to  obtain  0.1X  statistical 
repeatability  levels  with  1  to  5  seconds  per  angle  Increment  with  only  a 
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small  mass  density  of  added  seeding  scatterers  of  this  size  along  the 
measurement  path.  It  should  be  borne  In  mind  that  the  practical  limits 
of  flow  angle  precision  and  turbulence  Intensity  measurements  may  be 
set  by  data  rates  which  are  marginal  in  the  TGF  without  a  very  small 
amount  of  seeding.  The  achievement  of  optimal  seeding  will  require 
further  analysis,  design  and  experimentation. 

Measurements  in  the  Self-Adaptive  Wall  Tunnel  Indicate  that 
there  are  an  adequate  number  of  detectable  particles  In  the  facility 
for  rapid  LTA  measurements  (20-50/second)  at  Mach  0.9.  However,  due 
to  the  instability  of  the  particle  concentration,  higher  thresholds 
(with  lower  data  rates)  and  conditional  statistics  must  be  uBed  to 
reliably  measure  flow  angle  with  the  LTA.  Thia  can  probably  be 
accomplished  without  seeding.  However,  since  time  constraints  are 
more  severe  In  the  SAW,  it  would  be  time  efficient  to  seed  and  get 
the  measurements  down  to  one  second  per  angle.  This  would  allow  an 
entire  angle  scan  to  be  conducted  in  less  than  10  seconds.  If  this 
capability  were  coupled  with  a  faBt  computer  driven  traverse  mechanism, 
it  should  be  possible  to  obtain  very  precise  mean  flow  angle,  and 
turbulence  Intensity  at  5  per  run  minute.  Key  points  to  consider 
is  that  these  measurements  could  be  made  with  0.5  pm  diameter  particles 
or  smaller  which  are  responsive  to  flow  behind  shocks. 

The  overall  conclusion  from  this  series  of  measurements  Is 
that  the  LTA  system  is  capable  of  making  useful  aerodynamic  measure¬ 
ments  in  the  TGF  and  SAW  wind  tunnels  without  seeding.  These 
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capabilities  can  be  enhanced  by  seeding  with  very  minute  mass  quanti¬ 
ties  of  seeding  material.  Furthermore,  a  useful  data  base  was  pro¬ 
duced  in  the  TGF  with  the  microfighter  for  use  in  characterization 
of  the  effectiveness  of  the  adaptive  wall  concept. 
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LASER  TRANSIT  ANEMOMETER  WITH  MICROCOMPUTER 
AND  SPECIAL  DIGITAL  ELECTRONICS  I 
MEASUREMENTS  IN  SUPERSONIC  PLOWS 


by  W.  T.  Mayo,  Jr.,  A.  8.  Sure  And  T.  E,  Hunt 


Sp«ctron  Development  Laboratoriaa,  Inc. 

3303  Harbor  Boulevard,  Suita  G-3 
Coata  Maaa,  California  92626 

Summery  laprovement  >  to  both  the  optical  and  algnal  condi¬ 

tioning  ayataaa.  The  data  management  aye tan  (DNS) 
An  advancad  laaer  tranait  anemometer  ayataa  haa  bean  improved  under  SDL  private  funding.  Com- 

haa  been  developed  for  difficult  optical  flov  mat-  aarcial  eopiea  of  thia  even  acre  advanced  "fourth 

eurement  applicatlona.  Thia  co-axial  backacattar  generation"  ayatam  have  been  conatructed  and 

eyeten  maaauraa  scan  velocity,  flow  angle,  and  tur-  delivered, 
bulence  lntanalcy.  The  uaa  of  highly  focuaad  spots, 
apacial  pulee  detection  circuitry,  a  faat  aingle 
bit  digital  correlator,  and  a  microcomputer  data 
management  ayatam  allow  unaaedad  aubaonlc  meeaura- 
manta  to  rangaa  graater  chan  two  metera  and  aupar- 
aonic  measurements  with  naturally  occurring  aubmi- 
cron  partidaa  at  ahorter  rangea,  The  ayatam  haa 
excellent  wall  and  window  flare  rejection  capabili- 
tlea  and  may  be  uaed  for  boundary  layer  maaaura- 
aanta  near  aoltd  uodala.  It  la  compact  and  rugged 
and  haa  baan  damonatrated  for  auperaonlc  meaaura- 
manta  In  high  nolaa  and  vibration  environment*  at 
both  USAF  Arnold  Canter  and  the  NASA  Aama  Raaaarch 
Center.  Thia  LTA  ayatam  haa  alao  baen  damonatrated 
for  aubaonlc  axial  compraaaor  aaaaurananta.  In  re¬ 
cent  development*,  the  apaad,  aanaltlvlty  and  reli¬ 
ability  have  been  further  improved  in  preparation 
for  aaaaurananta  up  to  Mach  5  at  the  Arnold  Canter 
Von  Karaan  Facility  in  lata  turner  1979. 


1.  INTRODUCTION 


Laaer  tranait  aneaometry  (LTA)  or  "two-apot" 
techniques  have  recalved  attention  recently  due  to 
aanaltlvlty,  high  velocity  capabilities,  and  super¬ 
ior  flare  rejection  advent age* l-7.  An  LTA  ayataa 
aeeeuree  the  distribution  of  the  ties*  of  flight  of 
scattering  particle*  as  they  pa**  successively 
through  two  highly  focuaad  parallel  laaur  beaaa  aa 
illustrated  in  Figures  1  and  2.  The  appropriate¬ 
ness  of  using  pulse  canter  aatlaatlon  filters  with 
real-time  single-bit  croaa  correlation  end  the  need 
for  correlation  processors  with  S  to  10  nano-second 
time  resolution  for  supersonic  flows  have  both  baan 
docuaantad  recently®**. 

An  advanced  co-axial  backacattar  LTA  system 
wee  developed  and  daaonetrated  by  Bpectron  Develop¬ 
ment  Laboratories  under  eponeorehip  by  the  U.8.A. F. 
Arnold  Engineering  Development  Canter®,  NASA  Amt*®, 
end  private  funding.  Thle  eye tarn  included  a  third 
generation  corset  and  rugged  optica  head  with  a 
133  *M  laser,  optica,  spot-pair  rotation,  raal-tim* 
graphic  display  of  corralograma,  preliminary  on¬ 
line  date  analyeie  of  mean  epaed,  flow  angle,  tur¬ 
bulence  intensity,  end  meat  data  tranafar  or  stor¬ 
age  for  further  off-line  computer  processing. 

Mora  recently,  additional  development  funding 
from  the  U.S.A.P.  Arnold  Canter®  has  reaulted  In 


Pig.  1  -  Laaur  Beam  Geometry 

I  I 


Pig.  2  -  Operating  Principle  v*  •  d/u 
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Following  sect  ions  describe  various  aspects  of 
the  system  and  our  experience  to  date.  These  in¬ 
clude  the  system  configuration,  signal  conditioning 
electronics,  the  microcomputer  based  control  and 
data  acquisition  aubayatem,  and  dieplay  and  data 
reduction  algorithms.  This  is  followed  by  a  brief 
description  of  supersonic  measurements  and  conclu- 
eions  about  the  applicability  of  this  technology. 


During  our  LTA  development,  we  have  used  e  Mal¬ 
vern  K7023  30  ns  correlator  to  produce  a  single-bit 
cross  correlation  of  the  outputs  of  pulse-center 
estimation  filters  (discriminator*),  A  faster  "mul¬ 
tiple  etetlon"  correlation  proceesor  1*  presently 
under  development  *t  SDL  for  LTA  end  other  applies-, 
tion*.  Figure  3  is  an  example  crose  correlogrem  of 
transit  delay  intervals  from  a  AS  store  correlator. 
The  rme  turbulence  inteneity  for  this  example  is  7X. 


2.  SYSTEM  OVERVIEW 


Figure  3  is  *  photograph  of  ths  LTA  antics! 
head  and  tha  data  management  svatam  (DMS) .  A  digi¬ 
tal  cross  correlation  processor  (not  shown)  and 
interconnection  cables  complete  this  stand-alone 
system,  although  traverse-mechanism  and  oehar  inter¬ 
faces  are  available.  The  optical  head  contains 
aavaral  seta  of  printed  circuit  card  electronics. 
These  subsystems  include  a  rotator  prism  control 
card,  two  eats  of  filter  discriminators,  a  high  vol¬ 
tage  relay  circuit,  and  photomultiplier  tube  hous¬ 
ing  circuits.  Figure  A  la  a  photograph  of  one  of 
the  discriminator  cards. 


Fig.  I  -  LTA  Data  Management  System  and  Optical  Head 


Fig.  A  -  Example  of  Optical  Head  Electronics •' 
Discriminator  Card 


Fig.  3  -  An  Example  Cross  Corralogram 
with  IX  RMS  Turbulence  Intensity 


Tha  DMS  automatically  rotates  tha  two  focused 
basins  about  a  comnon  canter  through  a  aaquanca  of 
angles  with  tha  angular  step  else  determined  on 
line  from  the  turbulence  levsl  estimated  from  tha 
first  data  run.  Tha  correlator  is  started  and  the 
data  transferred  to  tha  DMS  in  proper  aequence. 

The  operator  ie  shown  a  graphic  display  of  tha  cor- 
re la grams ,  a  plot  of  signal  visibility  ratio  from 
which  mean  flow  angle  ie  extracted,  and  a  printout 
of  ths  estimates  of  mean  speed,  mean  angle,  and  ap¬ 
parent  turbulence  intensity. 


3.  OPTICAL  HEAD  ELECTRONICS 


Through  previous  experience,  on*  of  the  author* 
(Smart)  found  that  lsrga  electromagnetic  interfer¬ 
ence  fields  were  often  present  in  coann  test  envi¬ 
ronments  In  addition  to  high  levels  of  noise  and 
vibration.  Thus,  in  addition  to  designing  for  me¬ 
chanical  rigidity,  we  chose  to  locate  tha  analog 
electronics  inside  the  optic*  heed  where  short  low- 
impedance  leads  end  ihislJlng  could  be  used  to  stain- 
tain  speed,  sensitivity  end  good  lmunity  to  inter¬ 
ference.  Only  TTL  levsl  signals  are  transmitted 
over  th*  long  signal  and  control  cables  (30  meters 
or  acre) . 

Tha  photomultiplier  tube  circuit*  offered  quite 
A  challenge.  We  desired  good  single-photon  respone* 
(pulses  up  to  two  volts  in  msgnitudt  and  lee*  than 
10  ns  wide  with  little  or  no  ringing  and  e»th  good 
single  photoelectron  puls*  height  resolution)  while 
allowing  such  longer  classical  signals.  This  must 
be  accomplished  in  th*  pretence  of  well  flare  end 
bled*  flash  which,  for  brief  periuds,  drive  the  tube 
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output  far  bayond  the  allowed  average  /mode  current. 
The  PMTs  selected  were  12  stage,  2  Inch  diameter 
>311  9817B  linaar  last  tecuaed  tubes  with  S-20  pho¬ 
tocathode  material.  These  tubas  have  high  cathode 
maximum  current  limits  such  chat  the  cathode  is  pro¬ 
tected  by  the  Inclusion  of  a  series  100  MO  resistor. 
Good  single  photoelectron  pulse  tieight  response  and 
first  stage  gain  is  aasurad  regardlass  of  the  high- 
voltage  setting  by  two  series  150V  Zener  diodes  be¬ 
tween  the  cathode  and  the  first  dynoda. 

The  Inputs  and  outputs  arc  illustrated  in  Fig¬ 
ure  6.  The  gain  auppresrion  dooa  not  disconnect 


will  come  on  when  either  tube  is  tripped.  The  red 
light  is  incorporated  in  a  switch  which  resets  the 
relays  when  pushed. 

Thera  is  a  separata  discriminator  card  for  each 
of  tha  two  aaparata  photomultiplier  tubea.  Each 
discriminator  card  la  a  "mother"  board  with  six  sep¬ 
arata  plug-in  "daughter"  modules  as  illustrated  in 
Figure  4.  Each  module  ia  for  a  pulse  width  with  a 
dynamic  range  of  4tl.  Tha  ranges  overlap,  aa  they 
increment  in  powers  of  two  with  range  1  being  25  ns 
to  100  ns,  range  2  being  50  ns  to  200  ns,  etc.  Fig¬ 
ure  6  illustrates  one  of  the  plug-in  modules  funct¬ 
ionally.  The  lov-pass  filters  are  designed  to  have 


LIGHT 


CATH00E  PROTECT 
RESISTOR 


Fig.  6  -  Photodetection  Electronics. 


the  high  voltage  but  removes  tha  lnteratage  gain  at 
the  first  dynoda  thus  preventing  predictable  tran¬ 
sient  optical  overload  from  fatiguing  the  dynoda 
electron  emitteru.  Tha  PKT  anode  load  la  50  (1. 

There  are  two  direct  coupled  stages  of  preamplifi¬ 
cation  which  together  give  a  gain  of  x  100  and  do 
not  change  tha  polarity  of  the  output  algnal.  These 
preamplifiers  saturate  at  about  -  6  volts.  In  gen¬ 
eral,  the  uee  of  euch  preamplifiers  is  not  desira¬ 
ble,  and  for  wind  tunnel  applications  away  from 
model  surfaces,  a  14  stage  tube  with  one  lOx  preamp 
is  a  better  selection.  However,  the  applications 
of  interest  often  require  proximity  to  walla;  thus, 
the  preamplifiers  are  required  to  avoid  dynoda  fa¬ 
tigue  through  excessive  average  currant.  Each  pre¬ 
amplifier  output  drives  a  50  Q  co-axial  cable  termi¬ 
nated  on  a  discriminator  board.  Figure  7  is  a  pho¬ 
tograph  of  two  phototube  assemblies  which  are  par¬ 
tially  disassembled  to  show  the  circuitry. 

Between  the  two  preamplifiers  there  is  an  in¬ 
tegrating  monitor,  located  physically  on  tha  high 
voltage  relay  board,  which  detects  a  condition 
which  could,  if  sustained,  damage  the  later  dynodes 
or  anode  of  the  photomultiplier  tuba.  This  trip 
condition  corresponds  to  a  dc  erode  condition  of 
200  uA  for  a  period  of  0. 5s.  The  high  voltage  re¬ 
lay  responds  to  the  integrating  monitor  (actually 
an  active  low-pase  filter)  by  disconnecting  the 
high  voltage  i  o  protect  the  anode  and  later  stages 
of  tha  dynrfde  chain.  Bach  tube  has  a  separate  sense 
line  and  a  separate  high  voltage  relay  but  a  red 
overload  light  on  both  tha  optical  head  and  the  DM8 


Fig,  7  -  Photomultiplier  Tuba  Assemblies 


FB0H  PHT  ASSEMSLV 


Fig.  8  -  Schematic  of  Single  Filter 
Discriminator  Plug-in  Module. 
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nearly  syonatrlc,  non- ringing,  CauasLan  Impulse  ra- 
aponaas  whose  deration  ara  kT  and  4kT  (or  tha  de¬ 
tection  filter  and  the  pulae  ?«ntar  estimation  fil¬ 
ter,  raupactively.  Tha  rat.,*  aalactad  la  k,  from  1 
to  6.  Tha  1/s*  pulae  width  T  (or  tha  (aateat  range 
la'  23  na  (or  tha  ataadard  unite;  special  hlgher- 
apaad  plug-  in  modules  have  been  developed®  (or  ap¬ 
plication  up  to  Mach  10. 

In  operation,  the  uaar  electronically  aalacta 
tha  dlaci:ltainator  range  (or  which  aoat  o f  the  apot 
transit  durationa  T.  will  lie  within  the  range 
kt  <Tj<4kt  .  Under  these  conditions,  the  detection 
filter  (kx°)  has  a  wider  bendwidth  than  tha  daaai- 
cal  signals  to  be  detected,  so  there  is  no  statis¬ 
tical  biaaalng  due  to  attenuation  of  tha  (aster  pul¬ 
ses.  Once  the  detection  threshold  is  exceeded,  tha 
pulse-center  (liters  ara  enabled.  Since  tha  Inpul at 
reaponsa  o(  these  filters  is  longer  in  duration  than 
the  spot  transit  tines,  the  aingle-photon  pulaaa  ara 
smoothed,  even  whan  only  one  or  a  few  such  pulaaa 
occur  in  a  spot  transit  group.  The  coeparator  thus 
locates  tha  weighted  average  canter  of  a  photon 
pulse  cluster. 

The  discriminator  cards  include  a  second  level 
detection  circuit  just  following  the  buffer  anpli- 
flar  which  la  not  shown  in  the  simple  block  diagram 
This  circuit  inhibits  the  discriminator  outputs  if 
the  signals  excead  a  preset  level.  Tha  preset  level 
is  typically  -2.5  volts  to  avoid  very  large  inpulsas 
(frost  large  scetterers  or  other  causes).  By  adjust¬ 
ing  tha  PKT  high  voltage  aupply,  the  average  pulae 
level  out  of  tha  praanpliflers  la  eat  between  -0.2 
and  -2.0  volts. 

Tha  notor  control  board  supplies  power  to  the 
dc  geared  motor  which  drives  the  prism  rotator  such 
aa  to  aero  the  difference  between  the  logic  word 
provided  by  tha  tsicroprocaaaor  end  the  analog-to- 
dlgltal  conversion  of  the  signal  from  the  potentio¬ 
meter  connected  without  backlash  to  tha  prism  rota¬ 
tor.  for  incremental  angle  steps  lass  than  2  de¬ 
grees,  0.1  degree  precision  la  obtained  in  less 
than  0.3s.  The  value  of  ths  instruction  dispatch 
word  is  displayed  at  tha  microprocasaor  as  is  the 
word  which  corresponds  to  the  angular  position  ul¬ 
timately  obtained.  Tha  prism  motor  may  also  ba 
driven  in  either  fast  or  alow  mode  in  either  direc¬ 
tion  from  the  manual  switches  on  the  DH8  consols. 

In  this  mode  the  readout  of  actual  position  to  the 
DMS  screen  ia  sustained  and  updated  many  times  per 
second. 


4.  DATA  MANAGEMENT  SYSTEM  (DMS) 


The  DMS  (see  Figure  3)  is  normally  situated 
remote  from  tha  optical  head  by  up  to  30  a  and  con¬ 
tains  the  microprocessor,  keyboard,  visual  display 
unit,  two  mini  floppy  disk  units,  operating  soft¬ 
ware  end  a  number  of  manual  controla.  Figure  9 
shows  tha  labeled  front  panel  layout  which  ldantl- 
flea  tha  items  shown  in  the  photograph  of  Figure  1. 
Inside,  there  are  low  voltage  power  supplies  and 
other  electronic  hardware  to  augment  ths  micropro¬ 
cessor.  These  ere  48k  bytes  of  BAM,  inter faces, 
and  the  window  end  delay  unit  (cne  channel  of  which 
controls  tha  "live  tlms"  of  one  discriminator  and 
tha  other  ths  "blanking  time"  of  both  photomulti¬ 
pliers  for  periodic  flow  applications.  )  The 
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parallel  ports  (Interfaces)  provided  In  the  DMS  In¬ 
clude  a  buffered  8  bit  duplex  computer  Interface 
with  four  additional  control  Unas  and  a  48  bit 
correlator  lntarfscs.  Ths  correlator  ports  are  pro¬ 
grammable  with  raspact  to  tha  direction  of  informa¬ 
tion  flow.  Thus,  interfaces  to  a  variety  of  other 
computers  and  correlators  may  be  configured  by  soft¬ 
ware  without  hsrdvars  changes.  Sines  all  software 
is  provided  on  disk  rather  thsu  permanent  read-only 
memory,  ths  system  is  quite  flexible. 

Software  ss  provided  is  read  from  Disk  1  on 
startup  and  Disk  2  is  reserved  for  date  storage. 
Either  will  serve  both  functions,  if  necessary.  All 
Initiation  and  data  heading  functions  are  entered 
from  ths  keyboard  in  response  to  questlone  posed  on 
ths  screen  ss  illustrated  in  ths  example  given  in 
Table  1.  These  keyboard  entries  are  written  at  ini¬ 
tialisation  heading  before  each  data  record.  Fur¬ 
ther  details  concerning  ths  data  processing  algori¬ 
thms  are  presented  in  Section  3  below. 

Ths  photomultipliers  are  each  supplied  from 
chair  own  high  voltage  supply  so  that  fins  tuning 
on  gain  may  bs  performed.  The  remote  laser  start/ 
safe  la  s  removable  key  switch  with  laser-reedy 
indication.  This  prevents  ths  need  for  access  to 
tha  laser  power  supply  except  when  ons  wishes  to 
change  power  level e  or  the  control  mode.  There  le 
a  leser  output  potter  meter  on  the  control  coneole. 
Manual  controls  which  may  override  the  computer  ara 
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Table  1  -  Example  of  Dec*  Entry  Prior  to  Run 


QMS  RR0HF1  TO  user 


example 


•PUI  IN  OATA  DISKETTE*? 


Confirm  tint  this  Is  Oont 


b)  Give  »n  approximate  value  for  veloc¬ 
ity,  turbulence  and  angle  In  minimum 
time, 

c)  Print  these  estimates  on  paper  for 
brief  assessment, 


DATE  (WITH,  DAT) 

NO.  OF  STONES  •  ? 

SPOT  SEPAR  IISS)  •  ? 

NO.  OF  It  ADCS  •  T 

HINDOO  WIDTH  FOR  DATE  SIGNAL  •  7 

DEIAT  FOR  DATE  SIGNAC  •  ? 

WINDOW  WIDTH  FOR  ILANK1N3  SIGNAL  ■  ? 

OCLAT  FOR  DLANKINQ  SIGNAL  ■  1 

NO.  OF  POINTS  TO  USE  ON  EACH  SIDE 
OF  THE  PEAK  *  ? 

NO.  OF  OtCTA-T'S  (POWER  OF  10)  •  ? 


12,  25 

40  (40  Store  corrt lltor ) 

0.27 

29 

Lldth  snd  Delay  In  Nlcrottcondi 


7  (Sm  Sictlon  5  ) 
6 


X,  Y,  2  •  7 


10.25.  3.45,  5.00 


TIME  (HOUR.  MINUTE )  ■  7 

OELTA-T  (MICROSEC)  •  7 

0ES1RED  RANGE  ■  7 

THETA  •  [VALUE  DISPLAYED 
REAI  TIME  UPDATE 

DESIRED  THETA  ■  7 

THEtA  *  180.5 
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provided  to  make  setting  up  the  system  easier. 

There  is  a  switch  which  allows  slther  manual  or  cota- 
putar  discriminator  range  selection.  Mode  switches 
for  mstor  control  to  facilitate  alignaant  of  the 
spot  orientation  with  the  flow  are  provided.  Even 
under  manual  control,  the  correct  prism  orientation 
is  read  to  the  screen  of  tha  DM8.  The  threshold 
control  tor  ths  discriminators  seta  over  a  wide 
range,  the  lovel  which  changes  tha  ranges  ths  ratio 
of  Slagle  photon  events  to  multiphoton  particle 
events. 


5.  SOFTWARE 


The  optical  and  electronic  system  for  transit 
anaRK>mstry  basically  acquires  two  channels  o t  digi¬ 
tal  events  whose  only  information  is  carried  at 
event  times.  To  measure  the  distribution  of  event 
times,  it  in  optimal  to  cross  correlate  the  two 
channels  where  the  delayed  channel  corresponds  to 
that  whirh  a  particle  crosses  time.  This  subject 
has  bean  previously  discussed'.  Tha  subject  of  op¬ 
timum  dots  extraction  froe  corralograms  can  gst  vary 
detailed.  Here  we  will  explain  whet  we  have  elected 
to  uv  with  the  eyatern  to  data.  Furthar  aoftware  re¬ 
visions  ere  possible  either  by  the  user  or  froe  SSL 
on  s  new  minidisk.  Ths  present  software  functional 

a)  Maks  a  vlaibla  record  of  data  quality 
so  that  the  user  may  make  on-line 
Improvements  to  operating  procedure, 
threshold  end  voltage  settings, 


d)  Store  or  transmit  and  identify  all 
acquired  date  for  later  processing 
with  superior  algorithms  in  a  large 
computer. 

Figure  1  is  a  schematic  view  of  tha  sampling 
volume  and  shows  how  the  velocity  le  resolved  in 
tha  axial  direction  but  only  gives  rise  to  s  mea¬ 
surement  of  true  spaed  if  It  is  in  ths  pls.te  of  ths 
two  parallel  cylinders.  Mors  simply,  this  can  ba 
viewed  as  in  Figure  2.  Bacauaa  there  is  no  way  of 
identifying  a  specific  particle,  we  have  chosen  to 
dacoda  tha  data  on  a  statistical  basis.  By  comput¬ 
ing  tha  times  batwtien  all  possible  'A  than  S'  pairs 
such  that  tha  daisy  does  not  sxcaad  ths  displayed 
delay  time  of  tha  correlator,  than  it  is  possible 
to  obtain  a  peak  from  those  particles  which  crossed 
A  and  than  B  and  a  background  from  all  others.  Over 
short  delay  intervals  it  is  permissible  to  assume 
that  uncorrelatad  events  from  particles  which  cross 
only  one  spot  or  from  other  random  events  ere  Pois¬ 
son  distributed.  This  it  not,  In  general,  true  for 
longer  correlation  times;  and,  indeed,  provides  a 
method  of  assessing  turbulence  teals  by  correlating 
tha  existence  of  particle  concsntrstlon  fluctua¬ 
tions9.  Over  ths  short  tlma  involved  in  a  typical 
transit,  these  effuets  are  not  aumlfaat  and  it  is 
permissible  to  assume  that  ths  "background"  Is  flat. 
Figure  5  shoved  a  typical  corralogram  for  about  71 
turbulence  and  a  48  store  correlator.  Upon  this 
plot  era  marked  tha  various  operating  parameters 
which  are  used  by  ths  DMS  to  evaluate  tha  flow 
parameters. 

There  ere  many  simplifications  in  ths  follow¬ 
ing  analysis.  It  is  not  to  be  scorned  for  that,  as 
it  is  Intended  as  a  quick-look  device  and,  if  used 
correctly  with  proper  precautions,  is  cnpabls  of 
much  better  than  IX  accuracy  in  velocity,  particu¬ 
larly  for  low  turbulence,  and  no  moro  than  a  few 
percent  error  on  turbulence  for  turbulence  up  to 
about  13X.  For  very  low  turbulence  mors  stores  and/ 
or  smaller  &T  corrslstor  resolution  may  be  required. 
The  low  turbulence  limit  aleo  depends  on  the  degree' 
of  deconvolution  of  tha  instruiMne  broadening  func¬ 
tion'  one  may  wish  to  attempt  in  poet-detection 
processing.  Ths  standard  algorithm?  are  prssantsd 
below. 

Based  on  expsrlsnes  and  "a  prlerl"  knowledge  of 
ths  flow,  tha  operator  select*  e  numbor  of  scores 
HP  which  he  aspects  to  bound  the  hlgh^e:  turbuienoa 
conditions  of  the  corralogram.  Tha  region  R  is  than 
tha  sat  of  stores  located  .within  NT  stores  either 
side  of  ths  peak  IKX  of  tha  corralogram.  The  re¬ 
gion  I  le  all  the  rest  of  the  store*  a*  shoYm  in 
Figure  5.  The  first  step  in  the  data  reduction  Is 
to  compute  the  arithmetic  average  of  the  values  *c- 
cummulatsd  in  ths  stores  of  region  I.  This  1*  an 
aeciMts  of  tha  "whits"  level  due  to  uncorrelated 
background  pulses.  This  value  is  subtracted  from 
the  corrslogram.  Thus,  If  Y(l)  are  ths  corralogram 
vtluaa,  ve  obtain  I 

(1) 
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where  Q  la  the  total  "area"  under  the  peak  (see  Fig¬ 
ure  10)  which  normalizes  the  function  aa  a  probabi¬ 
lity  density  function  estimate;  i.e., 

(3)  Q  -  £t(I)-B 
H 


Fig.  10  -  Definition  of  Vlalblllty  Ratio 


the  baseline  la  considerable  as  illustrated  In  Fig¬ 
ure  5,  then  the  turbulence  Intensity  estimate  be¬ 
comes  significantly  affected  by  the  noiae_on  tho 
baseline  and,  hsnca,  the  choice  of  R  and  R.  Under 
such  conditions,  mors  data  ahould  bs  acquired  or 
sera  elaborate  curve  fit  parameter  extraction  algo¬ 
rithms  should  ba  employed  for  optimum  ratrlavsl  of 
information  from  the  data. 

The  DM8  uses  the  aatlsMte  of  turbulence  inten¬ 
sity  to  compute  an  incremental  angle  AS  for  spot- 
pair  rotation  which  proceeds  automatically  unless 
another  A6  Is  chosen.  The  computation  ofA0  assumes 
that  ths  turbulence  Is  isotropic  to  obtain  a  crude 
estimate  of  a  "good"  A6.  Thus,  If  the  number  of 
angles  at  which  observations  are  to  bs  made  la  N, 
than  AO  la  computed  aa  In  Table  2. 

Table  2  -  Selection  Formula  for  AO 


N  DT 

3  -  1.5  x  TI  x  180/n 
5  -  1.0  x  Tl  x  180/n 
7  -  0.7  x  Tl  x  180/it 


The  first  order  mean  speed  estimate  (assuming  that 
the  spot  pair  axis  la  directed  somewhat  along  the 
mean  flow  direction)  la  obtained  from  tha  mean  de¬ 
lay  number  KBAR  quite  simply  as 


W 


V  « 


d 

(KBAR  •  At) 


Unless  Tl  •  0,  In  which  case 

3  -  1.5  x  0.02  x  180/n  »  1.72 

5  -  1.0  x  0.02  x  180/n  -  1.15 

7  -  0.7  X  0.02  X  180/n  -  0.80 

and  Tl  is  displayed  aa  £  0. 


where  d  la  tha  physical  dlatanca  between  the  spots 
and  At  Is  the  correlator  clock  Interval.  This  est¬ 
imate  of  maan  spaed  baaed  on  the  inverse  of  the  mean 
transit  time  is  similar  to  tha  maan  velocity  compon¬ 
ent  estimate  obtained  from  burst  counter  processing 
of  f rings  LV  data  without  inversion.  Aa  is  wall 
knowrA  ,  the  use  of  tha  inverse  of  tha  mean  transit 
time  provides  first-order  correction  of  the  particle 
rate  blessing.  Thus,  for  low  turbulanr.s  flows  tha 
mean  speed  estimate  is  excellent.  At  higher  turbu¬ 
lence  levels,  say  greater  than  101,  higher  order 
data  processing  with  tha  inclusion  of  the  "pit 
slice"  effect  (accepted  flow-angle  at  any  one  angle 
depends  on  particle  else  tor  a  fixed  threshold)  and 
other  effects  may  be  Included  In  a  more  precise 
model. 

For  simplicity  the  present  turbulence  inten¬ 
sity  estimate,  TI,  is  calculated  as  a  root  maan 
square  deviation  uf  the  transit  time  distribution 
instead  of  the  velocity  distribution.  This  may  bs 
improved  in  ths  future.  In  terms  of  the  effective 
width,  F,  of  the  corrslogram  peak  we  obtain 

(5)  TI  »  F/2KBAR 

where  this  somewhat  unusual  notation  ie  used  in  tha 
definition  of  the  signal  visibility  ratio  function 
below i 

(6)  F  -  2\/£  (1-KBAR) 1 

The  abova  formula  for  mean  speed  la  quits  in¬ 
sensitive  to  tha  valuaa  of  NF  chosen.  If  tha  data 
acquisition  time  is  so  short  that  ths  scattat  in 


During  tha  automatic  angle  incrementation  se¬ 
quence,  tha  signal  "visibility  ratio,"  VR(L) ,  is 
computed.  The  angle  la  0 i  +■  LAB  where  0j  is  the 
initial  angle  selected  manually  or  via  computer. 

The  quantity  VR(l.)  is  the  ratio  of  tha  arsa  under 
tha  delay  curvs,  Q,  and  the  sum  of  this  area  with 
that  portiou  of  the  background  within  the  "width" 

F  of  ths  correlogrami 

(7)  VRCL)  -  ^ 

These  quantities  are  illustrated  in  Figure  10.  Tha 
value  of  ths  function  VU(L)  lies  between  taro  and 
unity.  Tha  ssro  condition  is  obtained  whan  tha 
spot  pair  sxls  is  rotated  so  far  away  from  ths  maan 
flow  angle  that  no  corrslatsd  events  occur  and  Q 
goes  to  zero. 

Tha  mean  flow  angle  occurs  in  ths  vicinity  of 
the  peak  of  ths  plot  of  VR(L).  The  present  DM8 
software  uses  a  3-pol.nt  parabolic  fit  of  tha  peak 
and  tha  value  of  VR(L)  on  each  side  of  the  peak  to 
obtain  tha  Man  flow-angle  estlMte  as  illustrated 
in  Figure  11.  Ths  idea  behind  the  definition  of 
ths  visibility  ratio  function  is  that  it  removes 
tha  effects  of  particle  rats  fluctuations  when  ths 
threshold  is  sat  sufficiently 'high  for  all  detected 
pulses  to  be  from  particles.  At  lower  thresholds 
and  pa-ticularly  when  tha  probe  voIum  is  near  a 
flare  producing  surface,  many  of  the  detected  pul¬ 
ses  are  single  photoalectron  background  pulaes 
whose  rats  nay  be  subject  to  angular  orientation  of 
ths  spots.  Under  such  conditions,  the  raw  date 
rate  Q  which  has  background  subtracted  out,  My  bs 
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Fig,  11  -  Example  Screen  Plot  of  VR(L) 
and  Display  of  Reduced  Data. 

a  better  parameter  to  fit  for  mean  flow  angle.  All 
of  tha  flow  eetlmatlon  aoftware  la  written  In  BARIC 
to  allow  aaae  of  future  modification.  In  fact,  all 
of  tha  operating  software  la  provided  on  disk,  ra¬ 
ther  than  permanent  memory  and  can  easily  be  changed 


6.  MEASUREMENTS 


At  tha  time  of  this  writing  (July  197$)  addi¬ 
tional  supersonic  maaauramente  are  scheduled  to  ba 
conducted  at  the  Arnold  Engineering  Development  Cen¬ 
ter  Von  Karmaix  Facility,  but  these  will  not  ba  con¬ 
ducted  before  the  aubsiisaian  deadline.  These  taste 
will  be  through  windows  two  Inches  thick,  without 
seeding,  in  Tunnel  A  at  Mach  numbers  in  the  range 
1.3  to  5.0.  Tha  reeulte  of  these  teste  will  be  dee* 
cribed  at  the  Congraee  If  they  are  available  at  that 
time,  and  will  be  provided  later  In  other  publica¬ 
tions.  In  this  section  we  will  briefly  review  the 
supersonic  measurements  made  with  the  prototype  sys¬ 
tem  at  tha  Arnold  Centar,  September  197B,  and  at 
NASA  Ames,  October  1978.  Thase  experiments  are, des¬ 
cribed  in  more  detail  In  the  contract  reports5 •  . 

Tha  objective  of  the  AEDC  measurements  was  to 
demonstrate  an  LTA  system  with  a  digital  correlator 
and  to  compare  the  measurement*  with  elaultansouslj 
obtalnod  fringe  veloclmeter  measurements.  Ths  pro¬ 
totype  system  was  transported  to  the  Arnold  Center 
on  Suptxuiber  18,  1978.  Demonstration  aoaaursm'antt 
along  the  centerline  of  a  supersonic  Jet  with  a 
Mach  disk  were  nude,  with  the  LTA  system!  however, 
due  to  equipment  difficulties  experienced  by  the 
Arnold  Center  personnel  and  time  and  fundlnQ  con¬ 
straints  which  precluded  further  SDL  measurements, 
no  simultaneous  measuresaanta  ware  accomplished. 

The  AEDC  tests  were  conducted  using  a  free  Jet 
tast  (unseadad)  Installation  located  In  Propulsion 
Research  Cell  (R-1A-1)  of  the  Engine  Test  Facility 
at  the  Arnold  Engineering  Development  Canter.  As- 
aistsnee  wee  provided  to  SDL  personnel  by  Mr.  T.  V. 
del,  RTF,  Mr.  Virgil  Cline,  PVT,  and  the  facility 
technicians.  A  1"  diameter  underaxpandad  unhaetad 
Jat  with  a  praasure  ratio  of  A  was  used.  This  Jet 
has  been  tested  extensively  In  the  pest  with  both 
laser  valocimetry  end  shadowgraph  techniques  as  re¬ 
ported  by  Barnett  and  Giel1"  in  1977.  Tha  first 
and  second  'Mach  disks  were  located  1.73  Inches  and 
3.39  Inchss  downstream^®.  Due  to  other  activities 
which  had  occurred  since  the  calibrations  of  tha 


Jet,  a  right  anglo  bend  had  been  installed  in  the 
Jet  supply  pipe.  It  1b  thus  possible  that  differ¬ 
ences  in  the  flow  would  occur  with  respect  to  the 

earlier  measurements. 

The  output  lena  used  for  the  Jet  test  was  a 
672  ran  focal  length  lene  which  reeulted  In  a  throw, 
or  range,  from  tha  front  of  tha  inetrusmnt  of  597 
e  spot  separation  of  0.489  mm  as  calibrated  prior  to 
Inetrument  (hipping,  end  e  spot  else  of  approximately 
20  micrometera  diameter.  The  axial  speed  end  appar¬ 
ent  turbulence  intensities  as  determined  by  the  op¬ 
erator  end  the  date  management  syetem  ware  recorded 
for  several  locations  along  the  esntar  line  of  tha 
jat.  Flgura  12  Is  a  plot  of  tha  axial  apaad  mea¬ 
sured  along  tha  jat  cantarllna. 


Dovnitrue  Px«t«nc«  (Inchsi) 


Fig.  12  -  Axial  Speed  on  Jet  Centerline. 

Careful  study  of  all  the  AEDC  date  and  the  ap¬ 
paratus  Indicated  that  tha  selection  of  the  colli¬ 
mator  lens  for  ths  prototype  eyatem  wee  inadequate 
to  properly  focue  return  scattered  radiation  on  to 
the  fiber  optic  light  guides  employed  in  the  system, 
/he  resulting  loss  of  scattered  rsdistion  resulted 
in  s  receiver  efficiency  of  less  then  51.  Calcula¬ 
tions  had  indicated  that  measurements  from  parti¬ 
cles  down  to  0.3  micrometers  In  diameter  would  ba 
achievable.  The  particle  leg  evident  In  the  plot 
of  data  serose  tha  Mach  disk  indicated  Chet  the 
first  prototype  system  did  not  achieve  this  sensi¬ 
tivity.  Contractual  tima  and  funding  limitations 
prevented  further  data  collection  or  analysis  at 
that  time. 

Ths  purpose  of  ths  NASA  Ames  experiment  was  to 
obtain  data  for  design  of  s  larger  future  eyatem 
for  the  NASA  Aaxaa  11-foot  transonic  tunnel,  if  such 
a  system  appeared  feasible  from  tha  testa.  The  teat 
wee  more  difficult  then  ths  AEDC  teat  In  that  there 
wee  a  2-inch  thick  window  between  the  LTA  system 
end  the  flow  field,  end  tha  larger  pertlclea  impact 
out  of  a  recirculating  tunnel.  The  measurements 
ware  conducted  in  a  piggyback  test  al  the  6-foot 
transonic  tunnel  without  seeding  et  Mach  numbers 
from  0.6  to  1.6.  Repeatability  of  the  seen  flow 
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measurements  was  as  good  as  0.15X.  The  tunnel  waa 
not  seeded.  Minor  improvements  in  the  optical  sys¬ 
tem  had  been  made  to  improve  the  optical  efficiency 
of  tho  recelvere  to  nearly  10X  but  inaufflclant  time 
had  bean  available  to  replace  the  output  collimator. 
Computer  simulations  of  the  signals  using  Mia  scat¬ 
tering  codas  and  measured  optical  paramatara  indi¬ 
cated  that  data  was  obtained  from  particles  down  to 
about  0.5  micrometers  in  disaster,  but  no  standing 
shocks  wars  available  for  relaxation  tecta  in  this 
exercise. 

In  December  1978,  the  prototype  system  waa 
used  to  make  measurements  in  a  research  axial  coo- 
pressor  at  the  General  Electric  Co.  Soma  of  tha 
data  was  obtained  unaeedad  and  some  waa  obtained  us¬ 
ing  0.25  micrometers  diameter  solid  seed  material. 
For  that  teat,  which  will  be  reported  in  datail  at 
the  Spring  1980  A8ME  Conference  in  Naw  Orleans,  a 
diffraction  limited  collimator  Ians  waa  used  with 
marked  Increase  in  optical  efficiency  and  thus  sys¬ 
tem  sensitivity. 


7.  CONCLUDING  REMARKS 

An  advanced  microcomputer  automated  laser  tran¬ 
sit  anemometer  system  (LTA)  has  been  designed,  con¬ 
structed,  and  delivered.  This  LTA  system  has  been 
demonstrated  lo  ba  capable  of  coaxial  backacattar 
measurements  of  unseeded  flow  near  walla  and  sur¬ 
faces  and  at  velocities  up  to  Mach  1.6.  Further 
teats  at  apaade  up  to  Mach  5  are  planned.  Tha  pro¬ 
totype  system  has  been  transported  to  three  sepa¬ 
rate  real  world  tests  with  high  noiae  and  vibration 
environments  and  performed  well  in  all  cases.  The 
succeaaes  support  tha  design  choices  which  lncluda 
the  extensive  use  of  printed  circuit  electronics  in¬ 
side  tha  compact  optical  head.  Tha  usual  complaint 
concerning  data  collection  tine  with  an  LTA  system 
has  baan  substantially  overcome  by  microcomputer  au¬ 
tomation  of  the  syatam  control  and  tha  provision  of 
on-line  flow  field  parameter  aatlmation  and  graphic 
operator  oriented  screen  displays.  Tha  LTA  system 
is  more  aanaltive  than  frlnga  LV  systems  and  is  ex¬ 
pected  to  find  wide  application  to  measurement  prob¬ 
lems  requiring  good  sensitivity  to  the  cuhmicron 
particles  which  abound  in  moat  instances  without 
seeding.  A  practical  LTA  syatam  with  unique  fea¬ 
tures  has  been  developed  and  demonstrated, 
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Following  equation  (1),  please  insert: 


Next,  the  mean  delay  store  number  KBAR  is  computed 

(2)  KBAR  -  £  •  1 
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SPECTRON  LASER  TRANSIT  ANEMOMETER  DEMONSTRATIONS 
AT  ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
SEPTEMBER  1979 


•  UNSEEDED  FLOW 
ft  BACKSCATTER 

«  80  M  LASER  POWER 

•  %  MILLIMETER  RESOLUTION 


TUNNEL,  A  TUNNEL  B 

•  MACH  NO.  A, 5  ft  MACH  NO.  8 

.  725  M/S  •  1200  M/s 

•  5s  •  50  s 

•  0 Al  VELOCITY  PRECISION  •  0.6%  VELOCITY  PRECISION 

•  0.1*  ARC 
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A.E.D.C.  VON  KARMAN  FACILITY 
TUNNEL  A  SEPTEMBER  10,  1979 


DATA  COLLECTION  TIME 
MEASURED  VELOCITY 
CALCULATED  VELOCITY 
DIFFERENCE 
ANGULAR  DEVIATION 
EVENT  RATE 


5  S 

726,1  ms"-*-  ±  Q JfiZ.  r.m.s. 
720,7  ms-1 

o.m 

4  0.10°  ARC 
10/S 


SPECTRON  LTA 
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A.E.D.C,  VON  KARMAN  FACILITY 
TUNNEL  B  SEPTEMBER  1 1,  1979 


DATA  COLLECTION  TIME 

50s 

MEASURED  VELOCITY 

1205  MS~1  *  QiiS  R.M.Si 

CALCULATED  VELOCITY 

1184  ms“1  ±  QM.  RiM.Si 

DIFFERENCE 

1.77% 

EVENT  RATE 

10/s 

SPECTRON  LTA 
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TRANSONIC  COMPRESSOR 
NAVAL  POSTGRADUATE  SCHOOL  DEMONSTRATION 
SEPTEMBER  21,  1979 

NATURAL  PARTICLES  ONLY 
UNPREPARED  DIFFUSE  METAL  SURFACES 
5s  DATA  COLLECTION 
-150  EVENTS/MEASUREMENT 
SPEED  •  135, A  ms"1  ±  0.2% 

ANGLE  :  22,0°  arc  ±0.7°  arc 

TURBULENCE:  1.9  ±  0.2% 

DATA  RATE  =  30/s 

BACKGROUND  PHOTON  RATE  5  x  107/s 
SUBMICRON  PARTICLES 
FULL  BACKSCATTER 


SPECTRON  LTA 
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LASER  TRANSIT  ANEMOMETER 
DEMONSTRATIONS 
BACKSCATTER  UNSEEDED  IOOmW 


WHERE 

SEEED 

ms"1 

{-ARTICLES. 

DALLRATE 

~s 

(1) 

A.E.D.C.  1WJET 

A20 

LARGE 

2,000 

(2) 

N.A.S.A.  AMES 

6'  x  6# 

FOG/NATURAL 

10,000/200 

(3) 

G.E.  LiS.RiC. 

% 

wzu 

i 

100 

l 

(A) 

A.E.D.C.  VON  KARMAN 

Tunnel  A 

725 

CLEAN 

30 

Tunnel  B 

1200 

CLEAN 

10 

(5) 

Monterey) 

135 

NATURAL 

30 

SPECTRON  LTA 
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ACCOMPLISHMENTS 


•  ADVANCED  LTA  SYSTEM  DEVELOPED 

•  OPERATES  IN  DIFFICULT  ENVIRONMENTS 

•  MEASUREMENTS  IN  AXIAL  COMPRESSOR 

•  UNSEEDED  RACKSCATTER  UP  TO  MACH  8 
FROM  <UJ  VIA  PARTICLES 


SPECTRON  LTA 
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APPENDIX  B 


PAPERS  DESCRIBING  LTA  THEORY  AND  SYSTEMS 
A  TUTORIALS  LASER  TRANSIT  ANEMOMETRY 

A.  E.  Smart  and  W.  T.  Mayo,  Jr.,  "A  Tutorials  Laser  Transit  Anemometry," 
Proceedings  from  4th  International  Conference  on  Photon  Correlation  Tech¬ 
niques  in  Fluid  Mechanics,  Stanford  University,  Stanford,  CA,  25-27  August 
1980. 
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A  TUTORIAL:  LASER  TRANSIT  ANEMOMETRY 


A.  E.  Smart  and  W.  T.  Mayo,  Jr. 
Spectron  Development  Laboratories,  Inc. 
3303  Harbor  Boulevard,  Suite  G-3 
Costa  Mesa,  CA  92626 


1.0  INTRODUCTION 

In  the  early  part  of  the  last  decade  the  now  conventional  Fringe 
Laser  Doppler  techniques  were  gaining  acceptance  and  producing  useful 
measurements  in  a  variety  of  situations.  As  more  experience  became 
available  certain  areas,  for  example  turbomachinery  research,  proved 
to  be  difficult  applications.  After  a  brief  period  of  pessimism 
about  the  fringe  system  great  innovations  were  made  but  meanwhile 
other  techniques  were  sought.  Schodl  [15  through  18]  explored  and 
developed  an  idea  by  Tanner  [27,  28]  and  Thompson  [29]  which  changed 
the  illumination  in  the  sampling  volume  from  fringes  to  two  separated 
spots,  Figure  1.  Various  designations  which  are  all  synonyms  which 
reflect  the  choice  of  the  pioneering  researchers  include  laser  transit 
anemometry  (LTA) ,  transit  velocimetry,  (LTV),  "two-spot,"  and  "L2F." 


FIGURE  1.  ILLUMINATED  VOLUMES  a)  LDV,  b)  LTA 
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In  generalized  terms  transit  anemometry  (LTA)  may  be  regard  >d 
encoding  the  signal  differently  from  fringe  systems  (LDV)  «  .<3  vhe  two 
systems  may  be  regarded  as  Fourier  Transforms  of  each  other,  "ijvrc  2 
[22] .  Some  significant  differences  emerge  and  it  is  these  differences 
which  lead  to  preference  of  one  system  over  the  other  for  different 
applications.  Custom  suggests  more  differences  than  are  actually 
observed  when  the  parameters  are  tailored  to  a  specific  test  case 
but  in  broad  terms  the  significant  properties  are  summarized  here. 

The  comparisons  are  based  upon  a  theoretical  experiment  where  either 
fringe  or  transit  system  might  be  successfully  used. 

Details  of  optics  ara  given  in  Section  2,  information  retrieval 
is  discussed  in  Section  3  and  in  Section  4  we  return  to  the  quantifi¬ 
cation  of  those  factors  briefly  dibcussed  in  this  introduction.  An 
appendix  illustrates  some  of  i.hr  geometrical  properties  of  the  LTA 
sampling  principle. 

Salient  properties  which  are  in  3ome  way  special  to  the  LTA  rather 
than  previous  systems  are  now  discussed.  Increased  illumination  inten¬ 
sity  is  typically  obtained  with  the  LTA  optical  configuration  and  this 
leads  to  ability  to  observe  smaller  particlea  with  similar  laser  power. 
With  a  typical  cross  sectional  area  reduction  of  the  illuminating  beam 
of  100  to  500,  it  is  possible  to  restrict  the  sensitive  volume  length 
by  a  ratio  of  between  10  and  20.  This  reduces  the  distance  from  walls 
at  which  many  measurements  may  be  made  with  a  given  size  of  particle 
[14]. 


Because  the  volume  is  both  shorter  and  narrower  for  LTA  than  for 
comparable  LDV,  the  data  rate  may  be  leas  because  with  a  given  particle 
concentration,  size  distribution  and  flux  fewer  particles  give  rise  to 
a  detectable  signal.  This  is  partially  offset  by  the  natural  tendency 
for  increasing  particle  concentration  with  reducing  particle  size. 

With  the  LTA  system  the  signal  processing  by  photon  or  event  cor¬ 
relation  has  been  shown  to  be  vary  satisfactory  [2,  10,  11,  12,  13,  20, 
25,  26].  The  bunching  of  single  photons  from  particle  transits  is 
statistically  discriminable  above  the  background  single  photon  rate  by 
a  filter  and  thresholding,  resulting  in  separation  of  multiple  photon 
events  from  background  photon  flux  rates  greater  than  10.  The  prefer¬ 
ence  for  fast  digital  correlation  processing  is  the  reason  for  the 
inclusion  of  LTA  applications  in  this  conference. 

Precision  can  be  high  with  LTA,  particularly  of  mean  velocity  and 
flow  angle,  0.1%  and  0.1°  arc  are  directly  attainable  in  many  situationa 
with  few  special  precautions.  Turbulence  may  be  measured  in  a  different 
way  from  that  available  from  former  technologies.  The  measurement  of 
speed  and  ite  distribution  in  a  chosen  direction  gives  access  to  esti¬ 
mates  of  turbulence  in  different  directions,  even  to  the  assessment  of 
stresses.  Early  systems  tended  to  be  most  suitable  for  turbulence 
between  2%  and  15%  but  this  is  now  [14]  extended  down  to  a  small 
fraction  of  one  percent.  The  upper  limit  is  based  upon  long  data 
gathering  time  and  as  yet  inadequate  theory. 
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2.0  OPTICS 


The  optical  design  of  LTA  systems  call  for  more  care  than  that 
formerly  necessary  in  optical  anemometry.  The  sampling  volumes  are 
smaller  and  the  stray  light  stops  capable  of  better  rejection.  The 
optical  performance  must  approach  theoretical  limits  to  use  the 
potential  advantages.  Figure  3,  which  may  be  referred  to  throughout 
this  section  has  its  components  labelled  in  three  ways.  [A]  represents 
essential  components ,  [B ]  represents  highly  desirable  components  and 
[C]  represents  components  which  depend  on  the  application. 

2.1  Beamspllttlnq 

Several  types  of  amplitude  division  beamsplitter  (Figure  3)  are 
in  use,  the  two  most  suitable  for  LTA  being  the  Rochon  and  Wollaston 
polarizing  beamsplitters.  Both  units  produce  collimated  but  slightly 
diverging  beams  of  crossed  polarization  from  parallel  input  beam  but 
these  authors  see  the  synmetry  about  the  axis  of  the  diverging  beams 
as  a  slight  advantage  of  the  Wollaston.  Other  beamsplitters  are  in 
use  and  offer  advantages  of  adjustability  but  this  can  cause  mechanical 
inconvenience  or  unreliability.  The  collimated  beam  from  the  laser  is 
thus  divided  into  two  collimated  beams  which  may  be  brought  to  two 
focal  spots  in  the  plane  of  a  lens  corrected  for  infinite  conjugates. 
This  lens  essentially  performs  a  Fourier  transform  function,  the 
splitting  feoe  of  the  prism  being  the  same  optical  distance  behind 
the  back  principal  plana  of  the  lens  aa  1b  the  focus  of  the  two  spots 
in  front  of  the  tront  principal  plane.  In  fact  this  distance  may  be 
slightly  modified,  in  situations  where  ultimately  parallel  spots  are 
essential  in  the  test  space,  to  yield  spots  which  are  converging  at 
the  first  focus.  ThiB  allows  the  relay  lenses,  whose  longitudinal 
magnification  would  otherwise  create  diverging  spots,  to  image  the 
final  test  Bpots  as  parallel. 


The  imaged  spots  are  relayed  into  the  test  space  by  back-to-back 
collimator  lenses  which  are  free  from  spherical  aberration  over  the 
full  aperture.  Using  this  method  the  spot  size  is  determined  by  the 
diffraction  limited  image  of  the  GausBian  input  laser  beam  via  the 
magnification  derived  from  the  focal  lengthB  of  the  three  collimator 
lenses.  The  separation  is  controlled  by  the  divergence  of  the  splitter 
prism,  the  focal  length  of  the  focusing  lens  and  the  magnification  of 
the  two  output  lenaes.  The  spot  diameter  to  separation  ratio  is  thus 
easily  designed  to  be  a  chosen  value  but  its  change  is  not  trivial. 


Quant itativaly  if  the  spot  separation  iB  s,  the  beamsplitter 
divergence  is  B,  the  input  laser  beam  diameter  is  D,  the  half  width 
at  the  Gaussian  waist  is  r  focusing  lens  and  internal  and  external 
lenses  of  focal  lengths  f3,  fj»  «nd  f3  respectively  we  may  write 
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*  Vs 


(i) 
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Optical  Switch  Ic] 
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and 


r  => 


_2 

H 


(2) 


giving  a  spot  separation  to  diameter  ratio  of 

jL  ttBd 
2r  ’  4X 


(3) 


In  practice  small  imperfections  in  optical  components  may  conspire  to 
reduce  this  slightly. 

2.2  Spot  Rotation 

In  the  LTA  technology  it  is  essential  to  rotate  the  orientation 
of  the  two  spots  and  three  ways  of  doing  this  have  come  into  uset 
rotation  of  the  beamsplitter  and  the  stops  separately  but  coupled 
together,  use  of  an  image  rotator  prior  to  the  output  lens  and  rotation 
of  the  whole  machine. 

The  last  method  produces  excellent  mechanical  correspondence 
between  output  and  return  images  but  is  bulky  and  inelegant.  The 
first  two  lead  to  a  more  compact  design.  The  first  has  more  moving 
parts  with  opportunity  for  less  than  perfect  alignment.  The  second 
is  efficient  and  rather  expensive.  It  ia  usual  to  arrange  for  computer 
control  of  rotation  with  corrective  feedback  of  the  angular  position 
attained  which  may  be  to  a  precision  of  0.1 6  arc  without  great 
difficulty. 

The  important  criteria  associated  with  the  ability  to  rotate  the 
spots  are  preciaion/accuracy ,  exact  placement  of  stops  independent  of 
angular  station,  preservation  of  beam  equality,  freedom  from  aberrations 
and  stationarity  of  center  of  gravity  of  spots  (rotation  without  dis¬ 
placement)  .  With  the  first  system  described  above,  it  is  necessary  to 
compensate  the  laser  polarisation  so  that  the  beams  do  not  change  their 
relative  intensity  with  rotation.  In  the  second  case,  since  the 
Wollaston  prism  remains  stationary,  small  adjustments  in  its  angular 
orientation  may  be  made  to  render  the  spots  of  identical  intensity. 

Sinca  the  laeer  is  typically  plane  polarized,  a  45*  Inclination  of 
the  axis  to  this  plane  yields  50i50  partition  with  normal  polarizations. 
Method  three,  with  rotation  of  the  laser,  also  is  likely  to  cause 
trouble  with  larger  lasers  and  is  probably  a  simplification  of  the 
optics  at  the  expense  of  mechanical  difficulties  and  physical  size. 

A  solution  that  has  been  used  by  the  present  authors  is  a  modified 
Mirror  Dove  prism  [23]  which  has  an  f/4  or  better  aperture.  It  leads 
to  compact  overall  system  length,  introduces  no  aberrations  into  a 
diverging  beam,  is  optically  efficiant  and  is  insensitive  to  polarization. 
There  are  no  surfaces  normal  to  the  optical  and  mechanical  axis,  which 
is  good  from  the  point  of  view  of  introduced  flare,  and  there  are  no 
refractive  elements.  The  input  contains  no  critically  toleranced 
components  but  is  difficult  to  align  initially. 
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2 . 3  Stops 


Probably  the  most  important  design  property  of  LTA  systems  is  the 
precise  specification  of  stops.  Their  function  is  to  reject  light  from 
all  sources  except  that  scattered  from  moving  particles  in  the  control 
or  sampling  volume.  To  that  end  a  simple  pinhole  overlain  by  the 
conjugate  image  of  each  sampling  region  would  seem  adequate;  however, 
there  are  several  potential  sources  of  flare  to  be  considered.  Multiple 
reflections  from  within  the  optical  system,  even  with  well  coated  optics, 
are  a  source  of  scattered  light.  These  may  be  reduced  by  careful  consid¬ 
eration  of  all  components  and  by  'double  stopping1  where  the  image  stops 
are  themselves  reimaged  on  two  more  stops.  This  is  effective  in  reducing 
second  and  higher  order  scatter  sources,  and  has  advantages  not  only  in 
improving  the  rejection  of  flare  generated  from  within  the  optical  sys¬ 
tem  but  also  from  multiple  scatter  in  any  test  situation  and  beam  walk- 
off  sources  in  windows.  Typical  rejection  factors  so  far  achieved  have 
been  up  to  twelve  orders  of  magnitude  defined  as  photon  rate  from  input 
laser  to  detected  single  photon  rate  from  P.M.T.  It  is  not  necessary 
to  do  better  than  this  for  turbomachinery  applications  where  the  internal 
flare  is  dominated  by  the  scatter  from  walls.  (One  always  desires  to 
measure  closer  to  walls  until  the  flare  mechanism  overrides  the  useful 
data.)  It  1b  very  useful  to  improve  this  figure  for  wind  tunnel  measure¬ 
ments  where  there  is  no  significant  flare  generated  near  the  test  space, 
although  window  flare  may  still  limit  the  usefulness  of  the  so  called 
improvement. 

In  addition  to  the  field  stops  discussed  so  far,  two  other  types 
of  stops  are  necessary  to  the  system.  An  aperture  stop  is  automatically 
provided  by  the  outer  edge  of  the  final  lens  and  by  the  odge  of  the 
turning  mirror  which  separates  outgoing  from  incoming  light.  It  is 
necessary  to  extend  the  stopped  area  of  the  turning  mirror  to  prevent 
diffracted  flare  from  its  edges.  This  is  conveniently  done  by  two  more 
stops,  one  of  which  restricts  the  input  beam  to  an  area  slightly  smaller 
than  the  turning  mirror  and  one  of  which  occludes  an  area  slightly  larger 
than  the  turning  mirror  from  behind.  These  latter  stops  ere  not  purely 
aperture  functions. 

Yet  another  stop  is  essential  to  improve  wall  proximity.  Flare 
may  be  seen  by  one  photomultiplier  which  has  arisen  from  light  in  the 
other  channel  [26,  27],  This  may  be  prevented  by  stopping  part  of  the 
received  area.  There  is  a  reduction  in  optical  efficiency  of  the  machine 
and  this  stop  is  only  necossary  when  wall  proximity  is  important.  This 
is  executed  as  an  extension  to  the  stop  .immediately  adjacent  to  the  turning 
mirror.  Some  wall  proximity  tests  are  documented  in  Reference  14. 

2.4  Aberrations 


The  LTA  principle  relies  heavily  on  having  very  small  spots  and 
well  defined  images.  It  is  necessary  to  approach  diffraction  limited 
performance  as  closely  as  possible  and  in  this  respect  LTA  systems  are 
more  difficult  to  design  and  build  than  equivalent  LDV  systems.  Because 
this  system  has  a  small  field  only  spherical  aberration  is  of  consequence 
among  the  first  order  Seidel  sums.  It  is,  therefore,  of  extreme  importance 
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to  obtain  lenses  whose  design  conjugates  conform  to  the  overall  system. 
Infinite  conjugate  design  is  needed  for  the  three  essential  lenses. 
Chromatic  correction  is  not  needed  unless  the  system  is  to  be  used  with 
a  fluorescent  seed  medium;  however,  fluorescent  particles  tend  to  be 
of  the  order  of  1  ym  diameter  or  larger  to  scatter  a  significant  amount 
of  wavelength  shifted  light  and  will  not  be  small  enough  for  many  appli¬ 
cations.  The  peculiar  advantages  of  LTA  may  not  be  so  well  taken  if  it 
i 8  necessary  to  add  such  large  particles  to  the  test  flow. 

Projected  spot  quality  and  receiver  efficiency  are  both  irqpaired 
by  normal  and  more  by  oblique  incidence  through  a  thick  window.  This 
effect  is  studied  and  quantified  in  Reference . 26,  but  typically  up  to 
2"  thick  windows  with  less  than  10s  incidence  do  not  cause  serious 
deterioration  with  f/4  systems.  If  windows  up  to  4"  thick  are  essential, 
near  normal  incidence  is  desirable.  Curved  windows  are  not  advised 
unless  they  are  quite  thin,  or  of  large  radius  of  curvature  as  they, 
like  inclined  plane  windows,  may  introduce  unacceptable  tangential  or 
sagittal  astigmatism. 

2.5  Mechanioal 

In  a  system  which  relies  upon  optical  precision,  mechanioal  proper¬ 
ties  are  at  least  as  important.  Small  size  leading  to  rigidity,  and  low 
weight  leading  to  ease  of  handling  and  traversing  are  important.  Reten¬ 
tion  of  alignment  in  bad  thermal  and  vibration  environments  is  an  essen¬ 
tial  feature  of  optical  anemometers  in  many  applications.  It  is  both 
more  critical  and  more  difficult  to  achieve  with  LTA  systems  than  with 
other  configurations  because  of  the  relatively  tighter  tolerances. 

2.6  Photomultipliers 

The  requirements  for  photodetectors  [Figure  3]  are  complex  in 
that  it  is  desirable  to  have  single  photon  counting  resolution  with 
low  time  and  amplitude  fluctuation.  The  detector  should  not  be  damaged 
by,  nor  be  subject  to  a  significant  recovery  time  after,  overload.  The 
details  of  this  are  the  subject  of  paper  No.  10  in  this  volume.  It  has 
been  convenient  to  make  the  stops  to  be  the  ends  of  single  fiber  optics 
whose  flexibility  makes  small  adjustments  both  easy  and  precise  without 
moving  the  comparatively  bulky  photomultiplier  tube  assemblies. 

3 . 0  INFORMATION  RETRIEVAL 

The  optical  configuration  of  the  sampling  volume  for  an  LTA  system 
is  quite  different  from  its  counterpart  in  LDV.  While  it  is  obvious  that 
an  LTA  system  measures  t'  nsit  time  it  is  not  so  obvious  that  this  is 
also  what  is  measured,  tr  s  indirectly,  with  an  LDV  system.  Many  ways 
of  regarding  the  LDV  fringe  signal  have  been  discussed,  with  models  of 
varying  completeness,  but  a  time  measurement,  in  one  or  other  guise,  has 
always  been  available.  The  significant  differences  are  in  the  precision 
and  interpretation  of  this  time,  and  in  the  form  of  the  measurement. 

Again  trivially  the  LDV  measures  the  velocity  component  at  right  angles 
to  the  fringes  and  has  a  precision  controlled  by  the  total  amount  of 
light  and  more  de  jusly,  the  number  of  fringes.  This  is  for  a  perfect 
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(optimum)  signal  processing  technique.  The  issue  of  matched  filtering 
based  on  'a  priori'  data  is  thorny  and  much  argued  but  applies  also 
to  LTA  in  a  much  simpler  way.  The  precision  here  may  be  enhanced  by 
moving  the  spots  further  apart  but  this,  for  a  given  level  of  turbulence, 
restricts  the  number  of  events,  and  because  of  this  reduction  in  statis¬ 
tical  numbers  may  not  lead  to  a  final  improvement  in  overall  accuracy. 

Suffice  it  to  note  that  the  types  of  velocity  measurement  are  quite 
different  in  LDV  and  LTA  and  each  may  have  the  more  appropriate  output 
for  a  specific  application. 

3.1  Time  Detection  for  T '"A 

The  detection  of  the  best  possible  estimate  for  the  time  when  the 
particle  crossed  the  peak  (nearest  point  to  the  center  of  the  illuminated 
region)  is  derived  by  filtering  the  signal  with  a  Gaussian  filter  whose 
roll  off  is  well  shaped  and  sufficiently  close  to  the  estimated  velocity 
(within  a  factor  of  two  or  better)  to  induce  minimal  distortion  into  the 
assessment  of  peak  time.  Such  discrimination  is  described  elsewhere 
([12],  and  paper  10  of  this  volume).  The  peak  disparity  iB  a  function  of 
how  many  photons  are  received  and  detected.  For  one  photon  or  lesB  the 
spread  on  the  time  is  the  same  fraction  as  the  spatial  spread  on  the 
spot  (both  referred  to  the  time  or  space  separation  between  spots) .  As 
the  photon  numbers  increase  this  improves  dramatically  as  the  filter  haB 
an  optimum  tendency  to  take  the  'center  of  gravity’  of  the  photon  burst 
so  long  as  the  photodetections  occur  within  a  period  short  compared  with 
the  inverse  filter  bandwidth. 

For  a  few  photons  per  transit,  typical  discriminator  jitter  is  in 
the  region  of  1  ns  to  2  ns  for  a  total  transit  time  of  one  spot  of  15  ns 
and  is  proportionally  greater  for  longer  transit  times.  This  mechanism 
has  a  similarly  good  performance  whether  the  number  of  photons  is  just 
three  or  four  or  corresponds  to  the  quasi-classical  continuum  of  photons 
resulting  in  a  pulse  of  widely  variable  height  [11].  To  achieve  the 
best  overall  performance  of  such  a  discriminator,  it  is  desirable  to 
separate  the  event  detection  filter  from  the  pulse  center  estimation 
bv  about  a  factor  nf  four.  This  allows  a  smaller  bandwidth 
on  the  event  detection  filter  without  compromising  center  finding 
sensitivity  and  precision  with  a  large  dynamic  range  of  spot  transit 
times  [13]. 

Discriminators  such  as  outlined  need  only  produce  a  clean  pu3.se  at 
a  fixed  time  after  the  estimated  peak.  Any  delay  is  not  important  so 
long  as  it  is  the  same  for  the  complete  detector/discriminator  units 
for  each  spot.  A  good  way  to  confirm  this  is  to  make  a  velocity  measure¬ 
ment  from  spot  A-*-B  (say)  and  then  rotate  the  spots  and  measure  from  &+>». 
Now  the  two  channels  may  be  checked  out  'in  toto*  or  at  various  sub¬ 
sections  by  exchanging  connections  at  different  points  in  the  system. 

3.2  Instrumental  Broadening 

There  1b  a  group  of  parameters  which  may  affect  the  observations 
and  it  is  useful  to  summarize  these  and  consider  (without  proof  here) 
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the  consequences  on  a  qualitative  basis.  Considering  a  possible 
trajectory  across  the  region  illuminated  by  the  two  spots  there  are 
several  ways  the  assumed  geometry,  Figure  4  and  6  is  not  satisfied 
by  the  real  world.  The  parameters  may  be  summ.irized. 

1)  Particles  miss  the  center  of  the  first  illuminated 
region  in  either  of  two  directions,  along  or  normal 
to  the  optical  axis  by  distances  by  and  b^. 

2)  Particles  may  not  move  at  right  angles  to  the  optical 
axis  by  an  angular  deviation  $  nor  in  the  common 
plane  of  the  axis  of  the  two  spots  by  an  analogous 
angular  deviation  8. 

3)  Particles  may  be  of  different  sizeB  (radius  «  a)  . 

Within  limits  defined  by  the  spot  diameter  to  separation  ratio  and 
the  increase  in  this  ratio  with  progression  along  the  optical  axis  in 
either  direction,  the  values  of  the  above  parameters  may  take  a  wide 
range  of  values.  For  isotropic  particle  loading,  the  positional  coordi¬ 
nates  under  1)  will  be  uniformly  and  equally  probable.  For  the  angular 
parameters  under  2)  it  is  necessary  to  anticipate  the  turbulence  distri¬ 
bution  in  velocity  space  since  this  is  convolved  with  the  geometrical 
transfer  function  of  the  whole  illuminated  space.  The  consequences 
under  3)  are  determined  by  the  particle  size  distribution  which  is  most 
often  not  known. 

The  detailed  examination  of  these  factors  show  that  in  any  real 
flow  situation  the  mean  transit  time  is  obtained  from  all  those  particles 
which  fall  within  the  boundary  conditions  defined  by  the  above  variables. 
The  interrelationships  are  complex. 

If  the  signals  are  such  that  only  one  photon  arrives  from  the 
experimental  configuration  at  one  time  then  the  broadening  may  simply 
be  described  as  the  product  of  the  sensitivity  quadric  (the  illumination 
intensity  times  the  receiver  aperture  function)  and  the  particle  path 
as  a  sampling  function.  This  arises  because  the  probability  of  scatter¬ 
ing  (and  hence  receiving  and  detecting)  a  single  photon  on  a  statistical 
basis  is  predicted  by  the  classical  continuum  descriptors.  The  instru¬ 
mental  broadening  is  thus  described  by  the  defined  illumination/reception 
function  and  the  integrated  probability  of  finding  particles  with  each 
of  the  five  parameter  restrictions  just  described.  If  there  is  some 
method  of  looking  at  a  large  number  of  photons,  say  enough  to  approximate 
the  classical  estimation  under  all  conditions,  then  it  is  possible  to  use 
the  same  peak  time  detection  to  estimate  the  center  of  the  classical  pulse. 
The  positions  in  the  x  direction  of  these  centers  do  not  change  with 
changes  in  b|_,  b//  or  a.  Hence,  there  is  no  intrinsic  error  nor  broadening. 
With  respect  to  0,  any  variation  will  be  simply  resolved  into  the  measure¬ 
ment  direction.  For  most  commonly  used  configurations  this  is  true  pro¬ 
vided  that  0<  ^70°.  This  is  in  accord  with  the  transit  system  resolving 
velocity  components  in  the  plane  of  the  two  spots  into  a  direction  parallel 
to  their  common  normal. 
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FIGURE  5 


ILLUMINATION  AND  RECEIVING  VOLUME 
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The  final  parameter  0  has  complicated  functional  properties.  If 
the  spot  diameter  to  separation  ratio  is  small,  say  a  few  percent, 
then  to  first  order  this  is  not  significant  but  for  high  precision  it 
should  be  considered,  particularly  in  turbulent  flow  situations  where 
measurements  are  made  with  the  spots  at  different  rotations  compared 
with  the  mean  flow,  i.e.,  various  discrete  values  of  0.  A  complete 
analysis  of  the  effects  of  these  five  parameters  is  not  within  the 
scope  of  this  tutorial  but  a  brief  examination  of  the  effects  and 
properties  of  a  typical  system  is  presented  in  the  Appendix. 

The  spot  transit  broadening  may  be  removed  by  peak  detection  in 
the  direction  of  the  mean  velocity  [14].  Where  measurements  of  turbu¬ 
lence  are  to  be  made,  it  is  apparent  that  the  broadening  at  right 
angles  to  the  mean  flow  may  be  more  difficult  to  remove.  The  success 
of  its  post  observation  deconvolution  will  depend  on  the  spot  separa¬ 
tion  to  diameter  ratio  and  the  level  of  turbulence  of  interest.  For 
mean  velocity  and  longitudinal  turbulence  measurement  the  6  broadening 
is  not  of  primary  concern. 

With  the  availability  of  increased  computing  power  it  may  bs 
desirable  to  look  at  all  the  data  from  a  variety  of  angles  and  to 
reconstruct  the  event  distribution  in  real  apace  and  averaged  over 
time.  The  non-linear  transformation  into  velocity  space  may  be  performed 
without  too  much  theoretical  objection  although  the  old  controversy  of 
the  particle  vs,  fluid  velocity  field  occurs  once  more.  It  is  now  in 
a  slightly  different  guise  because  of  the  topics  mentioned  earlier  in 
this  section  and  we  have  every  confidence  that  papers  in  this  conference 
will  illuminate  the  subject  further. 

If  an  array  of  data  points  is  available  in  velocity  space  then  the 
properties  of  the  flow/particles  may  be  quantified.  Two  provisions  are 
that  each  angle  Bhould  be  measured  over  a  time  which  is  statistically 
stationary  to  some  prescribed  variance  of  the  flow,  and  that  it  should 
comprise  sufficient  particle  samples  to  maintain  a  chosen  and  accepted 
maximum  variance  at  the  limiting  area  of  interest  of  the  final  velocity- 
space  distribution. 

Because  there  is  no  way  that  event  counts  far  from  the  mean  values 
can  satisfy  the  criteria  postulated,  it  may  be  good  either  to  fit  a  flow 
model  with  'a  priori'  reality  conditions  or  to  truncate  the  data  field 
at  given  distance  from  the  areas  of  maximum  interest.  Failure  to  do 
this  must  lead  to  catastrophic  breakdown  of  any  computing  algorithm, 
analogous  to  attempting  super-resolution  in  a  space  where  the  normalizing 
eigenvalues  are  already  hidden  beneath  the  statistical  noise  level.  This 
means  that  with  a  given  number  of  observed  events  only  certain  questions 
may  be  asked  of  the  data. 

4.0  QUANTITATIVE  PROPERTIES 

It  is  instructive  to  quantify  the  general  remarks  of  the  introduction. 
Many  approximations  will  be  made  to  show  trends  with  which  more  specific 
parametric  analyses  may  be  made  by  the  exacting  user. 
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4.1  Illumination  Intensity 


Referring  to  the  illumination  intensity,  the  overall  profile  is 
Gaussian  at  right  angles  to  the  Poynting  vector.  For  a  fringe  system 
this  is  oversimplified  as  this  region  is  illuminated  by  two  nearly  plane 
wavefronts  inclined  at  a  slight  angle  which  specifies  the  'fringe'  spacing. 
For  comparisons  this  model  is  adequate.  The  volume  size  need  not  'per  se' 
be  controlled  by  the  mechanism  of  LDV  vs.  LTA  but  in  practice  the  volume 
for  LTA  in  high  speed  applications  is  much  smaller  than  that  for  LDV. 

This  constraint  is  applied  because  the  LDV  area  must  typically  contain 
several  fringes  and,  depending  upon  the  application,  these  must  be 
widely  enough  spaced  not  to  exceed  the  frequency  handling  capability 
of  the  signal  processor  for  a  typical  high  speed  flow.  Related  limita¬ 
tions  do  exist  for  LTA  but  the  spots  may  be  made  the  sise  of  a  typical 
single  fringe  width  and  the  resolution  is  typically  much  better.  If 
the  velocity  is  low  enough  and  a  large  angular  aperture  is  available, 
a  fringe  system  may  show  axial  light  rejection  comparable  with  the  LTA, 
because  the  signal  processing  frequency  constraint  is  easod  and  very 
small  probe  volumes  become  tolerabla.  In  such  applications  as  auto¬ 
motive  reciprocating  engines  with  high  turbulence  and  low  velocities, 
no  immediate  advantages  are  seen  from  the  LTA  system  properties. 

Illuminated  spots  may  be  10  ym  to  20  ym  with  sufficient,  closely  spaced 
fringes  to  enable  excellent  experiments  to  be  performed  in  high  swirl, 
high  turbulence,  non-stationary  enviromnants . 

In  environments  where  velocities  are  high  and  the  physical  dimen¬ 
sions  of  the  test  assembly  are  not  proportionally  increased,  the  high 
frequencies  demanded  by  the  fringe  system  ere  not  so  easily  managed, 
in  turbomachinery  for  example.  Here  there  are  advantages  in  using 
smaller  particles  which  will  more  readily  follow  flowB  and  may  naturally 
(or  artificially)  be  present  in  larger  quantities.  We  will  handle  the 
smaller  size  first.  The  photon  flux  through  the  sampling  volume  ist 


F  ■  - j  photons/e/  unit  area.  (4) 

hwr 

Where  P  is  the  laser  power,  r  is  the  radius  and  hv  is  the  energy  per 
optical  photon. 


The  photon  flux  scattered  by  a  particle  of  radius  a  during  its 
transit  of  the  peak  is  ,rra2Fjq(a),  whars  q  is  the  scattering  efficiency. 
Hence  the  total  number  of  photons  scattered  during  the  transit,  G^,  will 
be  the  flux  multiplied  by  the  transit  time,  2r/V,  hence 


G  -  sai-la 

s  hvr  V 


(5) 


This  is  necessarily  approximate  because  no  allowance  is  made  for  which 
part  of  the  illuminated  region  is  intercepted  by  the  particle,  nor  the 
profile  function,  nor  the  collection  and  detection  efficiency  for  scat¬ 
tered  photons.  It  is  assumed  that  these  are  at  least  closely  comparable 
for  LDV  and  LTA  eystema  and  will,  thus,  not  have  significance  in  their 
comparison. 
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For  a  particle  of  radius  a,  therefore,  the  number  of  photons 
scattered  during  a  transit  of  an  LTA  volume  compared  with  an  LDV 
volume  is  2rD/rT  (where  and  r T  are  the  beam  waist  radii  for  LDV 
and  LTA  respectively)  assuming  that  the  same  laser  power  as  in  the 
LDV  volume  is  equally  divided  between  the  two  LTA  spots 

The  peak  photon  rate  for  the  two  systems,  because  of  the  dif¬ 
ferent  transit  times  varies  as  Hrj/x^)^.  Thus,  both  photon  number 
and  peak  rate  are  increasing  from  LDV  to  LTA  if  rD  >  rT,  as  is  the 
usual  case.  If  a  given  number  of  photone  or  a  given  rate  exaeas  over 
background  is  ths  oritsrion  of  signal  detectability  than  the  seme  success 
will  bs  mat  by  smallsr  particles  in  LTA  than  LDV.  As  an  asida,  we  note 
that  the  quality  of  optias  necessary  for  typical  LTA  systems  is  signifi¬ 
cantly  higher  than  that  for  LDV  which  operates  successfully  far  away 
from  diffraction  limitad  conditions,  sxoept  in  situations  where  the 
flows  ere  very  alow,  and  high  spatial  resolution  is  demanded . 

I 

The  comparative  retioa  here  and  throughout  this  section  are 
summarised  in  Table  1. 


TABLE  1 

RATIOS  OF  LTA  VS.  LDV  FOR  VARIOUS  PROPERTIES 
OF  EQUIVALENT  SYSTEMS 


Number  of  photons  scattered 
in  transit 


Peak  photon  rate  detected 


Ratio  LTA/LDV 


Sample  Volume  length 

Data  rate  for  monodisperaion 

-4 

Data  rate  for  n(a)  *  a  (il  ,a,_  »  minimum 
detectable  radius) 


Data  rate  for  N{«)  (with  allowance  for 
q(a)  and  retrieval  efficiency) 


Typical  ratio  of  probe  volume  radii  for  high  speed  flow  rD/rT  is  10  to  20. 
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4 . 2  Sampling  Volume  Length 

The  sampling  volume  for  each  system  may  not  be  directly  compared  on  the 
basis  of  instrument  function  but  it  is  possible  to  compare  one  spot  of  the 
LTA  with  the  sampling  volume  of  the  LDV  for  considerations  of  axial  sensiti¬ 
vity.  This  parameter  is  not  quantifiable  with  any  great  precision  but  on  the 
basis  of  perfect  optics  it  is  possible  to  assume  that  axial  acceptanoe  length 
is  approximately  proportional  to  r.  For  proximity  to  walls,  all  other  things 
being  equal,  the  smaller  r  the  smaller  is  L  and  closer  to  walls  may  data 
be  obtained.  If  this  .-ere  the  only  factor  involved,  the  ratio  of  wall  dis¬ 
tance  for  LTA  over  LDV  w.vid  be  rip/rD.  However,  wall  proximity  includes 
rejection  of  background  flare  as  well  as  probe  volume  length,  and  this  is 
even  more  strongly  influenced  by  the  radius  r  and  the  field  stop. 

4.3  Data  Rate 

Data  rate  may  be  quantified  only  with  some  reservations.  In  a 
cloud  of  monodisperse  particles,  whose  size  is  sufficient  to  give  good 
measurements  with  either  system,  the  data  rate  is  proportional  to  the 
swept  volume.  This  is  the  cross  sectional  area  normal  to  the  flow 
multiplied  by  the  velocity,  i.e.,  2LxtJV,  where  N  is  the  number  of 
particles  par  unit  volume,  aaah  of  which  is  of .radius  a.  The  data 
rate  ratio  for  LTA  to  LDV  is  therefore  (r^/r^)  ,  again  because  L  *  r. 

In  practice,  natural  aerosol  distributions  exhibit  a  number  density 
distribution  approximately  proportional  to  the  inverse  fourth  power  of 
radius,  of  tho  form  of  Figure  6.  The  observation  that  LTA  will  see 
smaller  particles  than  LDV  for  comparable  laser  power  and  other  suitable 
conditions  means  that  more  events  will  be  observed.  Supposing  that  each 
system  will  see  up  to  and  including  the  maximum  size  present,  of  whiah 

those  particles  which  are  above  a  given 


a  4  da  (6) 


(7) 

The  ratio  of  data  rate  for  LTA  to  LDV  is  now  modified  to 
it^/cp)2  x  (aD/aT)3  where  aD  and  a^  are  the  smallest  particles  which 
may  be  observed  consistently  by  the  LDV  and  LTA  systems  here  postulated. 

Once  more  there  is  a  modification  to  this  because  as  the  particles 
are  smaller  so  also  their  scatter  cross  section  (scattering  efficiency) 
reduces  still  faster  (q(a)  on  Figure  6) .  If  a  given  number  of  photons 
is  necessary  for  detection,  we  will  assume  that  K  more  photons  are 
required  for  LTA  than  LLV  (Note  that  K  <  1)  then 


2 

Now  in  the  range  of  interest  0.1  |jm  to  1  Dm  the  value  of  q(«)  «  a  , 


there  is  a  negligible  number, 
size  a  will  be  given  by 

v*>  *  T 

N  (a)  «  a'3 
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FIGURE  6.  DEPENDENCE  OF  SCATTERING  EFFICIENCY  AND  CUMULATIVE 
NUMBER  OF  PARTICLES  OF  RADIUS  GREATER  THAN  a  ON 
PARTICLE  RADIUS,  a 
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Hence,  the  limiting  particle 
related  as 


sires  for  each  system  are  approximately 


M4  .  2_!i 

“D/  K  rD 


(9) 


which  defines  the  limiting  ratio  of  minimum  urfobla  particle  sizes  as 


*T 


(10) 


and  hence  the  equivalent  value  of  data  rate  ratio  becomes 


r  r 


This  is  only  true  if  the  velocity  direction  produces  no  angular 
excursions  outside  the  range  of  acceptability  of  the  two  spots  (low 
turbulence) .  The  ratio  reduaas  rapidly  when  turbulence  increases  such 
that  only  a  small  number  of  particles  produce  events  on  both  channels 
in  any  one  orientation.  The  above  ratio  tends,  therefore,  to  be  a 
'best'  achievable  with  natural  particulates. 

4.4  Detection 


Detection  is  another  complex  area.  In  principle,  the  detection 
of  a  transit  event  is  easier  than  the  Assessment  of  frequency  within 
a  longer  pulse,  but  as  each  system  approaches  the  regime  of  quantum 
realisation  then  other  factors  come  into  play.  A  typical  spectrum  of 
a  classical  fringe  crossing  with  added  noise  looks  like  Figure  7(a).  The 
component  near  zero  (the  intensity  fluctuation  spectrum)  is  usually 
rejected  by  a  band  pass  filter,  represented  by  the  dotted  window.  In 
an  ideal  fringe  processor  this  window  is  a  tracked  envelope  the  same 
width  as  the  spectrum  of  a  single  transit  and  infinitely  agile.  This 
is  seldom  achievable.  In  the  transit  system  of  Figure  7(b),  the  'matched' 
filtaring  is  now  only  a  low  pass  function  which  extends  slightly  beyond 
1 Aj  to  cover  fluctuations  in  transit  time.  In  principle  this  need  only 
include  the  equivalent  low  passed  Intensity  fluctuation  of  Figure  7(a) 
(because  the  single  spot  trensit  hes  no  modulation)  but  in  practice 
single  spots  are  typically  of  the  order  of  sise  of  s  single  fringe  and 
hence  cover  almost  the  asms  range  in  frequency  space  except  for  the 
low  cut.  Certainly  its  noise  acceptance  is  comparable  and,  hence, 
appears  to  have  no  immediate  advantage. 

There  ia  an  advantage,  however,  and  this  may  best  be  seen  in  the 
time  domain.  Although  thia  may  not  be  seen  from  Figure  7,  it  is  clear 
from  Figure  8  that  the  non-linear  thrasholding  offers  an  immediate 
mechanism  for  rejection  of  'noise'  while  observing  'events'  whera  the 
diffarsnes  is  defined  ea  transient  photon  rate.  The  explanation  of 
this  principle  was  presented  in  reference  13. 
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b)  LTA 


FIGURE  7 


SPECTRAL  DISTRIBUTION  OF  LDV 
AND  LTA  SIGNALS 


In  an  experiment  without  additional  background  light  where  only 
light  scattered  by  isolated  particles  is  detected  it  has  been  shown 
that  30  to  60  photons  are  necessary  f or  1%  velocity  resolution  iB  a 
fringe  system  [4].  This  becomes  2  to  4  per  spot  for  a  transit  system. 
Indeed,  it  may  be  only  1  per  spot  if  the  spots  are  far  enough  apart  and 
there  are  no  other  photon  sources.  The  chosen  properties  of  the  filter 
for  the  fringe  system  are  not  so  necessary  for  the  transit  system  because 
any  orror  which  may  accrue  from  imporfcctionB  of  the  filter  will  bo  the 
same  for  both  opots  and  not  incur  any  first  ordor  orror. 

The  idea  of  using  an  appropriate  low  pass  filter  is  very  efficacious 
in  noise  rejection  because  of  the  possibility  of  analog  thresholding  to 
decide  whether  an  event  hae  occurred  [See  Figure  8].  A  single  photon 
(noise)  ideally  will  not  rise  to  this  threshold  when  convolved  with  a 
Gauesian  analog  filter  but  two  or  more  photons  occurring  within  the 
inverse  filter  bandwidth  will  cross  this  level  [Figure  8b].  In  any 
real  transit  experiment  the  peak  photon  rate  is  almost  always  too  high 
for  quantum  resolution  and  tha  analog  filtsring  represents  an  excellent 
method  of  discriminating  'events'  from  'no  events'  where  an  event  is 
now  associated  with  a  group  of  photons  from  a  single  particle.  A 
separate  filter  which  then  makes  an  estimate  of  time  of  occurrence  of 
the  highest  instantaneous  photon  flux  [Figure  8c]  may  generate  a  digital 
signal  [Figure  8d]  to  be  correlated  with  a  similar  system  which  views 
another  spot.  The  actual  fixed  delay  manifest  between  Figure  8c  and 
Gd  is  not  shown. 

In  practice  most  transit  systems  rely  on  having  2  or  3  or  more 
photons  scattered  from  a  particle  before  that  crossing  is  deemed  e 
successful  event.  This  is  s  property  of  where  the  systems  are  used 
end  mostly  their  advantages  are  great  when  operating  close  to  walls 
or  other  flare  sources.  Xn  these  situations  there  iu  a  background 
rate  of  detected  photons  which  are  almoat  always  individually  realised 
in  a  system  of  sufficient  speed  to  look  at  transits  with  nigh  precision. 

We  will  quote  ar.  example. 

At  400  ms  1  crossinq  a  10  ym  diameter  spot  takes  25  ns.  We  would 
look  at  the  signal  with  a  low  pass  filter  whose  Gaussian  shape  and  cut¬ 
off  point  will  yield  an  analog  pulse  two  to  three  times  greater  for 
three  photona  during  the  transit  than  for  one  at  any  other  time.  In 
principle,  a  background  rate  where  no  more  than  two  photons  occur  in 
25  ns  will  yield  no  events.  This  is  not  true  in  practice  because  of 
the  Poisson  nature  of  arrivals,  the  pulse  height  variations  of  the 
photomultiplier  (Paper  No.  10  In  this  volume],  and  other  faotors  but 
typically  three  photon  events  may  be  discriminated  in  the  presence  of 
background  of  up  to  107  photon j/s  of  'noiae'.  This  property  also 
removes  any  need  for  low  noise  photomultipliers  but  does  demand  a  photon 
counting  tuba  with  pulse  height  stability.  Again  with  some  considerations 
ths  speed  of  the  photomultiplier  end  detection  filter  circuits  control 
the  speeds  measurable  by  a  given  transit  system.  Basically  the  circuitry 
must  be  fast  snough  to  accomodate  signals  whoss  speed  depends  on  velocity 
and  optical  geometrical  specifications.  The  authors  have  measured  Mach  8 
(1200  ms~l)  with  12  ym  spots  to  batter  than  2%  using  15  ns  filters  [25]. 
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d)  Output  Pulaaa 


FIGURE  8 


DISCRIMINATOR  OPERATIC 
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Where  there  may  be  no  significant  flare  sources,  as  in  wind  tunnels 
away  from  models,  the  threshold  may  be  lowered  to  include  every  photon 
but  this  then  increases  the  'instrumental  broadening1  to  a  fractional 
level  corresponding  to  the  spot  diameter  to  separation  ratio.  This 
broadening  is  progressively  reduced  aa  the  number  of  detected  photons/ 
transit  increases. 

4.5  Data  Processing 

Several  methods  of  measuring  event  times  are  in  use,  the  two  most 
significant,  since  it  is  intervals  which  contain  the  required  information, 
are  direct  timing  and  correlation.  In  sparse  particle  situations,  where 
the  transit  time  from  one  spot  to  the  other  is  short  compared  with  the 
mean  interval  between  particles,  there  appears  at  first  to  be  no  signifi¬ 
cant  advantage  of  one  method  over  the  other.  Where  particles  may  occur 
more  frequently,  the  efficiency  of  timing  falls  rapidly  compared  with 
correlation.  In  all  cases  it  is  shown  very  much  more  aatisfactc'y  to 
use  two  photomultipliers.  With  one  PMT ,  timing  can  only  stop  on  tne 
next  consecutive  event,  and  this  yialds  a  sloping  background  which 
must  be  compensated  and  is  difficult  to  model  theoretically.  Auto¬ 
correlation  yields  a  flat  and  subtractable  background  but  this  is  four 
times  as  great  as  if  two  P.M.T.s  and  cross-correlation  are  used.  With 
two  P.M.T.s  and  timing  the  ’start1  may  be  on  one  event  on  one  channel 
and  the  ’stop'  may  be  on  the  next  event  in  the  other  channel.  A  more 
satisfactory  scheme  is  to  display  each  of  all  'stops'  on  the  second  chan¬ 
nel  which  occur  during  the  display  time  window  after  an  initiating  event 
on  the  first  channel.  The  latter  method  may  be  less  efficient  again, 
but  at  least  produces  a  subtractable  background.  These  authors  have 
found  in  practice  that  the  cross-correlation  principle  is  capable  of 
much  more  satisfactory  background  removal  and,  of  much  more  significance, 
of  operation  with  a  lower  threshold  of  particle  transit  event  detection. 

This  last  point  may  be  ejtplained  further.  The  typical  setting  of 
the  threshold  which  detscts  'real  events'  and  rejects  'noise  events', 
(particle  transits  and  background  noise  fluctuation  peaks  which  look 
the  same)  is  determined  by  the  need  to  obtain  an  'event  peak'  '■bove 
a  'noise  background'.  In  any  case  where  a  wall  is  approached,  >Ke 
operating  threshold  is  always  adjusted  to  boost  the  event  rate,  until 
the  noise  becomes  detrimentally  intrusive.  This  simulates  a  high  event 
rate  and  favors  cross-correlation t  specifically  this  may  be  seen  as 
the  capability  of  the  method  of  correlation  detection  to  look  at 
smaller  particles,  of  which  there  are  usually  many  more,  than  with 
timing  where  the  presence  of  background,  or  spurious  events  is  less 
easy  to  quantify  and  correct. 

With  either  detection  system  for  constructing  a  distribution  function 
for  event  interval  times  the  precision  of  the  final  estimates  of  mean 
velocit-,  turbulence  and  other  parameters,  such  as  flow  direction  and 
shear  stresses,  from  measurement  at  several  angles,  is  determined  by 
the  total  number  of  correlated  events  in  the  signal  peak  and  its  sepa¬ 
ration  from  the  background  noise.  If  the  background  has  a  statistically 
flat  probability  then  the  variance  of  any  assigned  'bin'  away  from  the 
signal  peak  is  equal  to  the  square  root  of  the  mean  occupation  number  of 
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the  background.  This  must  be  considered  when  using  the  occupation 
numbers  in  the  signal  peak  after  subtraction  of  the  mean  background 
computed  over  many  ‘bins'.  Schodl  [18]  has  suggested  a  scheme  for 
improving  estimates  of  mean  properties  where  several  angles  are 
measured.  Whereas  a  fully  resolved  veloaity  probability  distribution 
is  necessary  for  shear  stresses  the  accumulation  of  separate  ensemble 
averages  in  time  and  angle  means  that  average  properties  only  may 
be  obtained  to  precision  equivalent  to  the  former  unaveraged  measure¬ 
ments  with  only  one  tenth  of  the  data.  This  shortens  the  equivalent 
observation  time,  and  is  a  much  more  efficient  way  of  using  the  data. 
The  method  is  directly  applicable,  however  the  distributions  are 
obtained,  by  timing  or  by  correlation.  Depending  on  the  way  that 
correlation  is  calculated  the  distinction  between  limiting  and 
correlation  become  a  little  smeared. 

4.6  Precision 


Precision  is  not  such  a  contentious  subject  as  accuracy,  but  any 
fringe  vs.  transit  discussion  is  open  to  argument  on  either  topic.  It 
has  been  ouv  experience  that  LTA  systems  are  potentially  superior  in 
both  areas  but  that,  as  always,  thore  is  a  price.  The  first  and  most 
obvious  price  is  the  need  to  rotate  the  spots  in  time  and  thus  to  lose 
time  resolution  and  to  presume  statistical  atatior.arity  of  the  flow. 

Much  has  yet  to  be  said.  So  far,  less  creative  problem  solving  has 
been  addressed  to  the  LTA  to  extend  its  powers  in  the  way  that  has 
been  done  over  many  years  with  LDV. 

The  only  real  calibrations  for  LTA  are  laser  wavelength  and  spot 
spacing.  The  former  is  known  adequately  and  the  spot  spacing  may  be 
measured  to  better  than  0.1%  with  nothing  more  elaborate  than  a  polarising 
filter  and  a  tape  rule.  Calibration  for  LDV  is  more  difficult.  As  an 
aside,  in  areas  where  the  beams  are  projected  through  a  turbulent  region 
the  fringe  contrast  deteriorates  very  rapidly  and  a  small  change  in 
spacing  may  occur,  in  the  LTA  the  spots  may  broaden  but  are  unlikely 
to  incur  systematic  error  because  the  illuminating  beams  essentially 
travel  through  the  same  paths  in  the  fluid  until  they  begin  to  separate 
close  to  the  sampling  volume . 

Precision  of  mean  velocity  may  be  made  arbitrarily  high  by  moving 
the  spots  further  apart  and  by  taking  a  longer  time  average.  In  the 
caee  of  extremely  low  turbulence  this  is  sspecially  effective.  If  the 
experiment  time  is  shortened,  fewer  particles  will  cross  both  spots  and 
give  rise  to  correlated  data.  From  direct  statistical  considerations 
and  the  finite  bin  width  of  correlators,  precision  is  impaired  when  only 
a  few  events  are  available. 

High  angular  precision,  which  for  LDV  requires  crossed  fringes  and 
has  time  rewolved  properties,  is  achieved  by  spot  rotation.  This  is  easily 
achieved  to  0.2*  arc  for  low  turbulence  flow  and  flow  direction  has  been 
measured  by  observing  event  number  statistics  as  the  angle  is  scanned  either 
side  of  the  direction  corresponding  to  the  highest  correlated  event  rate. 

A  peak  fit  will  yield  the  mode  direction.  Because  of  the  finite  spot 
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diameter  to  separation  ratio,  this  angular  probability  curve  appears 
convolved  with  the  instrument  broadening  function  associated  with  the 
direction  window  width  within  which  trajectories  may  still  be  correlated. 
This  instrumental  broadening  is  of  the  same  order  as  the  longitudinal 
effects  before  reduction,  due  to  peak  time  detection  and  other  correlator 
defined  algorithms. 

Using  all  the  correlated  data  at  a  range  of  angles  the  statistical 
properties  of  the  flow  field  may  be  accumulated  in  time  and  inverted 
into  velocity  space  to  fit  known  or  hypothetical  functions  for  access 
to  higher  moments  of  velocity  and  stresses.  There  is,  so  far,  little 
work  on  this  as  the  instrument  functions  based  on  optical ,  geometrical 
and  signal  retrieval  phenomenology  become  very  complex  once  the  basic 
simple  formalism  is  extended.  For  all  practical  purposes  of  mean  flow 
velocity  and  angle  those  sophistries  art  of  no  importance  for  low 
turbulence.  In  that  mode  also  the  turbulence  may  be  well  measured 
up  to  (say)  ten  percent.  Above  this  level  and  referring  to  higher 
moments  or  cross  statistics  the  modelling  is  as  yet  in  its  infancy. 

There  are  some  suggestions  of  systematic  disparities  from  the  simple 
Intuitive  model  at  increasing  turbulence  levels  and  we  look  forward 
with  interest  to  new  work  in  this  area. 

4.7  Turbulence 


Turbulence  has  been  touched  on  briefly  in  the  previous  section  and 
at  various  periods  LTA  has  been  regarded  as  suitable  or  not  suitable 
for  very  high  or  very  low  turbulence.  Currently  it  seams  unwise  with¬ 
out  exhaustive  justification  to  believe  LTA  data  above  15%  turbulence. 

It  is  certainly  possible,  and  fairly  practical,  to  acquire  data  at  levels 
much  greater  than  this  and  while  the  data  is  not  in  question  its  inter¬ 
pretation  most  assuredly  is.  To  use  the  simple  models  at  30%  turbulence 
is  most  unwise. 

Below  15%  down  to  about  2%  LTA  with  cross-correlation  seems  at  least 
as  good  as  other  optical  anemometry  techniques  and  depending  upon  the 
application  may  be  much  better.  Earlier  work  with  correlators  suggested 
that  the  finite  bin  width  was  no  impairment  to  ultimate  velocity  precision 
provided  enough  events  were  observed  to  satisfy  the  requirements  of 
conventional  statistical  estimators.  This  means  that  even  though  the 
correlator  time  bins  were  quite  large  compared  with  the  time  resolution 
required  the  event  assignment  to  these  bins  is  very  precise  and  if  enough 
events  ocour  in  each  bin  then  interpolation  may  be  made  to  a  high  resolution 
which  is  a  function  of  bin  width. 

Racently  [14]  it  has  bssn  shown  that  there  is  a  first  order  inversion 
algorithm  for  subtracting  the  equivalent  mean  square  turbulence  apparent 
from  "bin  width"  effects  to  yield  the  ’true  turbulence'.  The  magnitude 
of  the  effect  depends  on  the  position  of  the  velocity  peak  with  respect 
to  the  time  boundaries  of  correlator  bins.  Using  the  presently  implemented 
scheme  turbulence  levels  below  0.5%  may  be  measured  with  repeatability  and 
close  correspondence  with  hot  wire*. 
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The  inversion  to  velocity  of  what  are  really  time  of  flight  measure¬ 
ments  is  well  understood  in  its  simplest  form.  The  measurement  of  mean 
transit  time  may  be  simply  divided  into  the  spot  separation  to  yield 
velocity.  Two  'errors'  [9]  which  compensate  to  first  order  are  involved 
in  this.  With  increasing  turbulence  this  compensation  deteriorates  and 
more  rigorous  algorithms  must  be  sought.  There  are  effects  of  "particle 
bias"  [8]  coming  from  the  lack  of  well  defined  correspondence  between  the 
particles  and  the  fluid  velocity  field  but  those  are  probably  insignifi¬ 
cant  below  5t  turbulence. 

5.0  CONCLUSIONS 

This  tutorial  is  a  brief  introduction  to  LTA  principle,  and  shows 
that  it  is  a  significant  alternative  to  the  widely  accepted  LDV  technology. 
The  two  areas  most  clearly  identified  to  date  are  turbomachinery  and  wind 
tunnols  where  it  may  have  important  advantages.  Table  2  is  a  brief  resume 
of  the  areas  and  propsrties  where  LTA  and  eingle  and  two  channel  LDV  may 
be  compared.  It  is  based  on  opinion  and  experience  to  data,  and  while 
any  item  may  be  the  object  of  an  interesting  discussion  the  broad  conclu¬ 
sions  reflect  opinion  trends. 

There  is  a  great  deal  of  work  to  be  done  on  the  interpretation  of 
the  data  derivable  from  LTA  instruments  and  this  confsrence  will  hslp 
to  fill  some  preeent  gaps  in  understanding. 
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TABLE  2 


Comparison  of  LTA  with  one  and  two  channel  LDV  systems  on  a  ' figure 
of  merit'  basis. 


Rank  Order  of  System  Superiority 


Experiment  Property 
(All  Others  Being  Equal) 

LTA 

LDV 

Single  Channel 

LDV 

Two  Channel 

High  Speed 

1 

2- 

2- 

Precision/Accuracy 

1 

3 

2 

Angular  Resolution 

1 

3 

2 

Low  Turbulence 

1 

3 

2 

Small  Particles 

1 

2- 

2- 

High  Turbulence 

3 

2 

1 

Data  Rate 

3 

1 

2 

Close  to  Walls 

1 

2 

3 

Cost 

2- 

1 

2- 

Ease  of  Use 

1- 

1- 

2 

Available  Information 

Time  Resolved  Data 

2 

1- 

1» 

1- 

2 

1- 

Spatial  Resolution 

1 

2« 

2- 
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>1-1.1 ....  .»* 


NOMENCLATURE 


a 


-  Particle  radius 


*L 


b 


// 


c 


f 


1' 


h 


-  'Impact  parameter'  normal  to  optical  axis  (y) 

[Sec  Figure  4] 

-  'Impact  parameter*  parallel  to  optical  axis  (y) 

[See  Figure  4]  2 

-  A  proportionality  constant  ■  i^xj 

-  Focal  lengths  of  focusing  lens,  internal  collimator 
and  external  transceiver  lenses  respectively 

-  Planck's  constant 


Ma) 

n(a) 

P<  ) 
q(a) 


t 


u,  v,  w 
x,  y,  a 


D 

F 

F. 


I 

K 


L 

N_(a) 

T 


-  Scatter  cross  section  allowing  for  the  efficiency 

-  Number  density  of  partiales  of  sisa  a  per  unit 
volume 

-  Probability  distribution  function  (different 
according  to  parameter) 

-  Partiole  scattering  efficiency  for  sisa  a 

2 

-  Beam  1/e  radius  at  Gaussian  waist  (sampling  volume) 

-  Spot  separation 

-  Transit  time  of  a  fringe  volume  or  between  spots 

-  Velocity  components 

-  Cartesian  coordinates 

-  Diameter  of  input  laser  beam  at  beamsplittar  face 

-  Photon  flux  density 

-  Photon  flux  through  sampling  volume 

-  Photona  scattarad  during  a  singla  transit 

-  Illumination  intensity 

-  Ratio  of  number  of  photons  needed  for  an  LTA 
detection  compared  with  LDV 

-  Sampling  volume  length 

-  Cumulative  number  of  particles  above  a  given  sisa  a 
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-  Laser  power 


P 

5  -  Scattered  signal  level  at  various  places 
-  Transit  tine 

V  -  Velocity 

3  -  Divergence  angle  of  beamsplitter  prism 

c  -  Ratio  of  sampling  volume  length  to  diameter  of  spot 

X  -  Optical  wavelength 

v  -  Optical  frequency 

6  -  Incident  angle  to  oonraon  plane  of  spots 

$  -  Incident  angle  in  common  plane  of  spots 

T  -  Transit  time  of  a  spot  or  fringe 

<  >  -  Ensemble  Average 


-  LTA 

-  LDV 
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APPENDIX 


GEOMETRICAL  PROPERTIES 


He  use  the  defined  notation  end  that  on  Figures  4  and  5  in  the 
text  and  the  classical  modelling  of  smooth  light  intensity  receivable 
from  all  points  on  the  trajectory.  This  is  proportional  to  the  product 
of  the  illumination  and  receiver  geometries  and  it  is  convenient  to 
express  the  measured  transit  time  from  peak  detection  methods  of  a 
single  particle  of  radius  a  as 


T 

a 


a  cos  6 

|  V]  ooa  $ 


(1) 


where  s  is  the  spot  separation  and  V2  »  u2  +  v2  +  w2,  $  and  6  are 
defined  on  Figure  4.  Since  it  is  only  possible  to  make  asaertiona 
about  the  known  direction  u,  as  that  is  the  only  measurement  available 
(the  spot  orientation  takes  the  coordinate  system  with  it  for  this 
analysis)  the  value  of  transit  time  beoomee 


T  -  —  COS20 

A  U 

If  the  flow  were  perfectly  aligned  with  x,  u  would  be  what  was  measured. 
If  it  is  not  then  what  we  really  measure  is  a  component  of  V  according 
to  the  sampling  and  geometrical  principles  of  the  LTA.  This  appendix 
discusses  the  disparity. 

In  any  real  experiment  the  value  of  Ta  must  be  accumulated  over 
many  particles,  a  statistically  representative  saxrple .  Hence,  the  mean 
value  of  Ta  is  obtained  by  integrating  over  the  boundaries  of  bj.,  b//, 

6  and  a.  He  will  postulate  functions  p(  )  (which  will  be  different 
for  each  parameter  in  parenthesis)  and  state  the  integral 


For  isotropic  particle  loading  p (bj_)  •  p(b//)  «  1  and  may  be  omitted. 

The  above  integral  is  not  simple  nor  separable  because  the  boundary 
conditions  are  not  mutually  independent.  Some  idea  of  this  is  available 
by  looking  in  detail  at  the  definitions  of  sampling  volume. 

The  illuminated  volume  is  Gaussian  in  radial  coordinates  from  the 
optical  axis  and  falls  off  approximately  hyperbolically  with  axial 


i 


distance.  This  may  be  represented  as 


The  fixed  dependence  upon  z  is  determined  by  the  definition  of  diffraction 
limited  spots.  The  receiving  aperture  is  not  so  easy  to  define  but  from 
Figure  5  it  may  be  approximated  as  a  hollow  cone  whose  inner  and  outer 
bounds  are  roughly  paraboloidal  and  hyperboloidal .  The  ratio  of  axial 
length  to  spot  diameter  is  a  function  of  the  design  of  the  receiving 
aperture  and  it  shall  be  designated  e.  For  a  practical  system  e  might 
be  ^50.  Because  of  this  '  t.  is  not  significant  to  modify  expression  (4) 
unless  the  flow  is  very  badly  misaligned  or  excessively  turbulent  neither 
of  which  conditions  is  under  discussion  here. 

Expression  <*)  must  be  extended  to  include  the  seoond  spot  under 
the  transform  x  •+  x-., 

I  (x.y.s)  12]  -  I<x,y,z)  +  I(x-e,y,z)  (5) 

If  a  particle  goes  through  these  spots  it  will  give  rise  to  a  signal 
which  has  dependence  upon  bj_.  b  >, ,  0,  a,  u,  v,  and  w  and  with  suf¬ 

ficient  devotion  the  express  dependence  may  be  written  down. 

With  borne  tedious  algebra  and  muoh  patience  it  is  possible  to  show 
certain  features. 

a)  The  variation  of  is  not  significant  if  peak  detection 
is  used  and  results  in  a  broadening  corresponding  to  the 
inarease  in  spot  diameter  with  s  if  quantum  resolved 
signals  are  processed. 

b)  The  variation  of  b^  is  interactive  with  6  and  will  result 
in  some  events  being  missed  in  peak  detection  if  bj.  and  0 
both  vary,  aB  it  is  almost  certain  they  will.  The  consequence 
is  complex  and  diacussed  further  under  d) .  In  quantum 
detection  no  events  are  missed  and  the  thresholding  effect 
resulting  in  a  sharp  boundary  is  transformed  into  a  steady 
Gaussian  roll-off. 

c)  Variation  of  $  haa  no  error  effeote  so  long  as  $  <  -  70*. 

The  velocity  component  in  the  zx  plane  is  resolved  into  the 
x  direction.  No  broadening  occurs  for  peak  time  detection 
but  significant  spread  may  bo  apparent  for  quantum  resolution. 
This  spread  will  increase  with  increasing  |$|. 
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d)  The  consequences  of  6  variation  really  get  complicated  1 

because  the  usual  mode  of  machine  operation  is  to  rotate  1 

the  spot  in  the  direction  6  to  make  measurements.  The  I 

boundary  on  0  is  dependent  upon  and  a  and  for  the  i 

first  time  a  systematic  error  in  Ta  may  occur,  bi  is  1 

less  significant  of  itself  but  it  could  increase  th®  1 

uncertainty  of  boundaries  in  0  and  a  domains.  For  peak  1 

detection  the  boundary  condition  on  a  defines  a  spot  I 

maximum  radius  of  intersection  beyond  which  no  event  | 

will  be  detected.  This  is  a  smooth  roll-off  with  quantum  I 

resolution  but  ia  still  an  attenuation.  For  a  constant  1 

photon  flux  there  is  also  a  reduction  in  photon  rate  with  J 

increase  in  total  speed,  jv|,  yielding  a  velocity  bias  1 

effect  which  is  a  boundary  definition  rather  than  a  1 

ubiquitous  scaling  parameter,  although  it  may  later  1 

cause  this  effect  because  there  are  two  spots.  It  is  | 

further  apparent  that  any  increase  in  ) 0 |  in  either  I 

direction  from  aero  has  an  effect  of  modifying  the  1 

measured  transit  time  in  the  same  direction.  This  1 

yields  a  systematic  erroneous  decrease  in  apparent  fl 

transit  time  average  if  there  is  a  wide  spread  of  0  H 

and  an  error  of  similar  sign  but  different  Bias  if  I 

the  0  variationa  are  not  symmetric  about  the  spot  to  1 

spot  direction.  This  effect  can  not  easily  be  quanti-  1 

fled  (which  means  1  haven't  done  it)  but  the  effect  1 

may  be  minimised  by  reducing  the  spot  diameter  to  9 

separation  ratio.  For  all  practical  application  for  1 

low  turbulence  this  error  seema  vanishingly  small.  I 

e)  Effects  with  variation  of  a  do  not  yield  systematic  )  I 

errors.  If  an  exclusion  limit  is  placed  upon  the 

maximum  particle  sise  then  particles  above  this  size  ! 

may  only  contribute  if  b^,  bi  ,  $,  6  or  combinations 

ere  increased  accordingly,  Tne  spread  will  increase  ! 

and  there  will  be  a  bias  against  those  particles  which  ; 

went  through  the  center  of  one  or  both  apots.  Smaller 
particles  will  have  a  smaller  range  of  the  above  para¬ 
meters  and  will  contribute  less  to  broadening.  The  1 

boundary  of  particles  too  small  is  a  roll-off  without 
error . 

| 

f)  The  effect  of  changes  in  u  is  to  change  the  equivalent  I 

signal  level  by  the  velocity  inverse  end  can  be  inverted 

from  average  time  measurements  by  a  simple  algorithm 

(s/Ta  for  single  event,  s/<Ta>  to  first  order  only  for  ) 

stationary  averages  in  low  turublence) .  For  peak 

detection  there  is  no  pystematic  error  and  no  broadening 

if  the  discriminators  are  perfect.  This,  naturally,  is 

only  so  for  measurement  of  the  velocity  properties  of  the 

particle  field  end  the  usual  corractiona  are  necessary  for 

the  trane formation  to  fluid  field,  together  with  the  usual 
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In  suninary,  the  various  geometrical  boundaries  are  reviewed  and 
none  is  found  to  yield  systematic  error  except  the  6  variation  which 
has  a  tendency  to  bias  the  measurements  of  transit  time  to  lower 
valueu,  leading  to  falsely  high  estimates  of  mean  velocity.  The 
magnitude  of  this  effect  is  very  email.  For  spot  separation  to 
diameter  ratio  of  25 tl  it  typically  might  be  of  the  order  of  1 
part  in  10*  or  lean  for  turbulence  sufficiently  high  to  offer 
transits  ovar  all  possible  6  values  and  could  possible  inarease 
to  3  parts  in  10*  for  turbulenca  of  less  than  0.2%  misaligned 
with  the  sensitive  direction  by  up  to  1*  (but  thought  to  be  cor¬ 
rectly  aligned) .  Any  angular  search  would  show  the  disparity 
immediately  es  it  beoomes  most  highly  sensitive  at  low  turbulence. 


The  systematic  errors  of  the  transit  geometry  are,  therefore, 

Bhovn  to  be  almost  insignificantly  small.  Any  residual  error  associated 
with  the  relation  of  the  particle  to  fluid  flow  field  must  be  considered 
separately.  Ha  ahow  our  final  Figure  A1  which  is  the  relationship 
of  the  boundary  contours  in  bj_  and  0  for  varying  values  of  the  parameter 
a.  This  is  no  more  than  the  flat  perspective  projection  of  a  three 
dimensional  section  of  the  8-space  of  coordinates  bj[,  b//,  <f>,  8,  a,  u, 
v,  w.  The  salient  features  are  briefly  as  follows. 

The  coordinate  boundaries  of  b /f  and  $,  together  with  v  do  not 
change  significantly  the  topographical  form  of  this  plot  nor  indeed 
do  they  produce  any  non-linearitiaa  which  would  render  it  lees  repre¬ 
sentative.  It  is  tschnioally  degsnsrste  under  w  as  no  measurement 
may  be  made  if  there  is  only  a  aignificant  w  component  (not  accommodated 
under  8).  The  v  component  is  similar  in  its  consequences  to  0,  but  with 
a  scaling  factor  which  would  not  show  on  this  plot.  The  affect  of  the 
magnitude  of  u  ia  to  scale  the  diagram  and  also  to  distort  the  a  coordi¬ 
nate  -  because  faster  particles  scatter  less  photons. 

The  remaining  three  parameters  are  plotted  as  boundary  values. 

Inside  the  shape  data  are  obtained.  Outside  they  are  not.  The  edges 
should  really  ba  represented  as  furry.  The  6  axis  is  shown  approximately 
quantitatively  with  an  axis  interval  of  0.5s  arc  for  a  25 tl  ratio  system. 
The  quantitative  marks  an  bj_  are  in  integral  multiples  of  rT  the  sensitive 
volume  minimum  radius  to  l/e?.  The  a  axis  is  not  typically  quantified 
but  scales  with  each  size  of  particle  and  with  inverse  of  speed  of  that 
particle.  The  quarter  space  thus  plotted  shows  simply  that  values  of 
b[,  6,  greater  than  certain  values  do  not  contribute  to  the  transit 
time  measurements  and  that  valuea  of  a  below  a  certain  sire  (or  for 
a  fixed  sire,  above  a  certain  speed)  may  not  contribute.  If  there  is 
a  method  of  excluding  particles  that  are  above  that  given  sire,  then 
the  represented  volume  develops  a  hollow  exclusion  core,  otherwise  the 
center  regime  also  contributes.  (leading  to  possible  slippage  errors 
from  large  particles  if  they  are  present.)  This  figure  must  be  convolved 
with  the  probabilities  of  finding  any  of  the  parameters  within  the  ranges 
of  interest  and  may  not  ba  shown  here.  It  would  correspond  to  a  distortion 
only.  For  example,  a  monodispersion  of  a  at  a  fixed  velocity  would  be 
a  ali.co  parallel  to  the  b^,  8  plane.  A  further  restriction  of  angle  8 
would  reduce  this  plane  to  a  horizontal  band. 
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